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Abstract: The effect of gravitational tidal forces on renormalized quantum fields 
■ propagating in curved spacetime is investigated and a generalisation of the optical the- 

^ ■ orem to curved spacetime is proved. In the case of QED, the interaction of tidal forces 

with the vacuum polarization cloud of virtual e^e~ pairs dressing the renormalized 
photon has been shown to produce several novel phenomena. In particular, the photon 
field amplitude can locally increase as well as decrease, corresponding to a negative 
imaginary part of the refractive index, in apparent violation of unitarity and the op- 
^ . tical theorem. Below threshold decays into e^e~ pairs may also occur. In this paper, 

these issues are studied from the point of view of a non-equilibrium initial-value prob- 
lem, with the field evolution from an initial null surface being calculated for physically 
distinct initial conditions and for both scalar field theories and QED. It is shown how a 
generalised version of the optical theorem, valid in curved spacetime, allows a local in- 
crease in amplitude while maintaining consistency with unitarity. The picture emerges 
of the field being dressed and undressed as it propagates through curved spacetime, 
with the local gravitational tidal forces determining the degree of dressing and hence 
the amplitude of the renormalized quantum field. These effects are illustrated with 
many examples, including a description of the undressing of a photon in the vicinity of 
a black hole singularity. 



1 Introduction 



Investigations of the effect of vacuum polarization on photon propagation in gravita- 
tional backgrounds have revealed many unexpected features of far wider importance 
for quantum field theory in curved spacetime. In particular, the geometry induces a 
much richer analytic structure for Green functions which, while preserving causality, 
fundamentally changes the assumptions behind established theorems in S-matrix the- 
ory and dispersion relations in flat space quantum field theory [1-5]. A variety of new 
phenomena associated with the lack of translation invariance also imply that many 
standard results based around unitarity and the optical theorem must be reassessed in 
curved spacetime. 

In QED, gravitational tidal forces act on the virtual cloud of electron-positron 
pairs which dress the photon with the result that photons do not simply propagate 
along classical null geodesies. Rather, the spacetime acts as an optical medium with 
a non-trivial refractive index n{u) [6-9]. Moreover, as originally found by Drummond 
and Hathrell [10], the low- frequency phase velocity arising from this effect may be 
super-luminal, which standard dispersion relations would imply is incompatible with 
causahty. 

In a series of papers [1-5], we have developed the theory of photon propagation in 
QED in curved spacetime and calculated the full frequency dependence of the refractive 
index in terms of a geometric quantity, the Van Vleck-Morette determinant, which char- 
acterises the null geodesic congruence around the classical trajectory. The important 
role of the Penrose limit [11-13] in identifying the salient features of the background 
geometry was identified and exploited to develop a rich phenomenology covering pho- 
ton propagation in a wide class of spacetimes. In particular, it was shown how the 
curved spacetime geometry induces a novel analytic structure for the Green functions 
and refractive index which modifies fundamental properties of QFT and S-matrix the- 
ory such as hermitian analyticity and the Kramers-Kronig dispersion relation [3]. This 
improved understanding of analyticity in curved spacetime allows a reconciliation of 
the apparent paradox between low-frequency super-luminal motion and causality. 

While this work resolved the problem of how causality is realised in QFT in curved 
spacetime, it led to a further apparent paradox, this time with unitarity. It was found 
that in certain backgrounds, including those associated with black holes, the imaginary 
part of the refractive index can be negative [3, 4]. This is in marked contrast to a con- 
ventional optical medium where Imn(a;) is always positive, corresponding to scattering 
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of photons from the beam and a reduction in the field amphtude. A negative lmn{u) 
would imply gain, with energy pumped in from an external source. In quantum field 
theory in flat space, the optical theorem relates Im?T,(ci;) to the rate of production of 
real e~^e~ pairs and therefore vanishes. Even if we considered an off-shell "photon" 
which is above the pair production threshold, lmn{uj) would still be manifestly posi- 
tive. In contrast, in curved spacetime we find examples where Im?T,(a;) can be negative, 
or non-vanishing and positive below the e~^e~ threshold. Clearly, understanding these 
phenomena and reconciling them with unitarity requires a careful reformulation of the 
optical theorem in curved spacetime. 

For instance, for QED in the weak curvature expansion we find that the refractive 
index matrix takes the form 



where u is the affine parameter along the photon ray and i,j = 1,2 are the two 
transverse space-like directions along which the polarization vector points. Riuju are 
components of the Riemann tensor and Ruu = R^mu- The first correction in (1.1) is 
the original Drummond-Hathrell result [10], which can lead to either a sub- or super- 
luminal low-frequency phase velocity. The second term, however, gives a non- vanishing 
contribution to lmn{u) and depending on whether the derivatives Riujuiu) are positive 
or negative along the photon trajectory can lead to an amplification or attenuation of 
the amplitude. 

In this paper, we investigate the realisation of unitarity in photon propagation in 
curved spacetime from the point-of-view of an initial- value problem. Light-front evolu- 
tion of renormalized quantum fields from an initial null surface is studied for a variety 
of field theories and initial conditions and the physical mechanisms responsible for the 
observed amplification or attenuation of the amplitude are explored. The centrepiece 
is the formulation of a generalisation of the optical theorem to curved spacetime. 

The key point is that in curved spacetime, the fundamental optical theorem relating 
a decay probability or cross-section and the imaginary part of a Green function must be 
viewed as a global result, integrated along the whole history of the decaying particles. 
The corresponding local identity, which is the standard theorem in fiat spacetime, 
loses the essential positivity needed for the identification with a probability. This 
vital difference between the local and global identities arises only in curved spacetime 
because of the lack, in general, of translation invariance along the particle trajectory. 
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Figure 1. Heuristic pictures whicli illustrate the behaviour of the renormalized photon 
made up of modes of the bare field (the wave) and virtual cloud of e~^e~ pairs. In the top 
diagram there are increasing Riemann tensor components (made precise in (1.1)) along the 
null coordinate u of the photon's propagation leading to an increase in virtual e~^e~ pairs and 
an attenuation of the photon modes (dressing). In the bottom diagram there are decreasing 
Riemann tensor components leading to the opposite effect and an amplification of the photon 
modes (undressing). 

As a result, the amplitude of a renormalized quantum field propagating through 
curved spacetime may be locally amplified, although integrated along its whole past 
trajectory there must be a net attenuation. This will be interpreted as a real-time 
dressing and undressing of the field by its cloud of virtual pairs. Locally the curvature 
can undress the field, resulting a reduced screening and an amplification of the renor- 
malized amplitude. Along its entire trajectory, however, the net dressing must remain 
positive. A heuristic picture of how dressing or undressing occurs is illustrated in Fig.l. 

These ideas are made precise in section 4, where an exact mathematical formulation 
of the optical theorem in curved spacetime is presented. The general theory is then 
illustrated with a detailed exploration of the initial value problem and field evolution in 
a variety of special cases, involving both different quantum field theories and different 
initial conditions, corresponding to a field in both the bare and dressed state. 
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The paper is presented as follows. We have already shown in previous work [3, 4] 
that for massless photon propagation through curved spacetime, the essential features 
of the curvature are captured by the Penrose plane- wave limit. We therefore restrict 
ourselves here to the study of quantum field theory in plane-wave backgrounds, al- 
though we shall consider both massless and massive theories. We therefore start in 
section 2 with a brief review of the nature of the Penrose limit and the geometry of 
plane waves. Section 3 sets up the initial- value problem and describes the different, 
physically distinct, initial conditions that we will consider. An essential feature is the 
Schwinger-Keldysh formalism [14-17] formalism for studying real-time, non-equilibrium 
phenomena in QFT and the relation of the causal field equations to the Feynman vac- 
uum polarization or self-energy is carefully explained. The general solution of the field 
equations for plane-wave backgrounds is then given and a re-summation technique, the 
dynamical renormalization group [18], is applied. The interpretation of attenuation 
and amplification of the amplitude as dressing and undressing of the renormalized field 
is introduced. 

The core of the paper is section 4, where the optical theorem in curved spacetime is 
presented. The relation of the local and integrated versions of the identity are explained 
in detail and the issue of unitarity and the positivity of imaginary parts of the Green 
functions and refractive index is addressed. The realisation of the optical theorem for 
different classes of initial conditions is also explained. 

The remainder of the paper illustrates these ideas for several quantum field theories 
- massless and massive scalar A(f)'^ theory in both d = 4 and 6 dimensions and QED. 
This allows us to address some special issues related to renormalization and also to 
discuss mass thresholds for decaying particles. 

First, we study real-time renormalization phenomena in flat spacetime, using a 
Laplace transform formalism, to develop insights into the initial-value problem and 
the role of initial conditions. This technique is then applied to QFT in symmetric 
(i.e. Cahen-Wallach [19]) plane wave spacetimes, which are translation invariant in 
the light-cone coordinate which serves as the affine parameter along the classical null 
geodesies. This reveals many interesting phenomena, including in some cases the below- 
threshold decay of the field A into (p pairs with a rate which is non-perturbative in the 
curvature. The consistency of all these results with our improved understanding of the 
optical theorem and the precise nature of the positivity constraint implied by unitarity 
is checked. 

Finally, we consider non-translationally-invariant backgrounds, in particular the 
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homogeneous plane waves which arise as the Penrose hmits of null geodesies in cosmo- 
logical or black hole spacetimes [4, 12]. In particular, we consider in detail the tidal 
effects on a renormalized photon field as it approaches a singularity and show explicitly 
how the photon can become undressed along such a trajectory. 

The insights and techniques developed in this paper should also have applications 
in a wider field. The methods of section 6 have already been used [18] to study real- 
time relaxation problems in quantum field theory, which are of relevance to inflationary 
and early-universe cosmology where the non-equilibrium evolution of scalar fields is im- 
portant (see e.g. [20-22]). Plane waves and Cahen-Wallach spaces are also important 
as backgrounds in string theory and have an important role in the AdS/CFT corre- 
spondence (see e.g. [23-25] and references therein). The study of QFT on plane waves 
has been pioneered in [26] and some more recent work motivated by string theory ap- 
pears in [27, 28]. Our work will answer some of the questions regarding the stability of 
interacting QFT in certain plane wave backgrounds posed in these latter two papers. 

2 Penrose Limit and QFT on Plane Wave Spacetimes 

To begin, we review some essential features of the geometry of plane wave spacetimes 
and how they arise in connection with photon propagation as Penrose limits. For a 
complete discussion, see our earlier papers, especially refs.[3, 4]. Some preliminary 
results on the formulation of QFT on a plane wave background are then described. 

2.1 Penrose limit and the geometry of plane wave spacetimes 

To illustrate how the Penrose limit arises and motivate our specialisation to plane wave 
spacetimes, consider first an interacting scalar field theory with an A(f)'^ interaction. 
The A field (playing the role of the photon in QED) is massless while the field (the 
electron) has mass m. Including the one-loop vacuum polarization, the equation of 
motion for the A field in the presence of a source J is 



We will define this precisely in the next section as an initial value problem, using the 
Schwinger-Keldysh formalism to identify n( X, X ^ clS the retarded vacuum polarization 
and A{x) as the 'Hn-in" VEV of the quantum field. 




(2.1) 
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We consider two limits: (i) u ^ and (ii) L<}\ ^ Ac, where u is the frequency 
of the wave solution for A, = ^ is the Compton wavelength of the "electron" , and 

is a typical curvature scale. The first is the geometric optics condition, which 
allows the description of the propagation of the free A field in terms of individual 
rays ("photon trajectories") following null geodesies. We then introduce Penrose, or 
adapted, coordinates (u, V,x"), a = 1,2, in which u is the affine parameter along 
the null geodesic, V is the associated null coordinate, and x° are transverse spatial 
coordinates. The free wave solution is then 

= (7(x)-^/'e'"''(") , (2.2) 

where g = —detg^i, and = d^'d{x) is the tangent vector to a null congruence 
g{k,k) = 0. Along the geodesies '(9(x) = is constant while the affine parameter 
u varies; V specifies the individual geodesic in the congruence. 

In our previous work [1-4], we have shown, using both worldline and conventional 
QFT methods, how the second limit allows us to simplify the evaluation of the vacuum 
polarization corrections to (2.2). To leading order in fH/m^ (where ~ is a typ- 
ical curvature scale), the vacuum polarization is determined by geodesic fiuctuations 
around the photon trajectory; in turn, this is determined entirely by geometric quan- 
tities, specifically the Van Vleck-Morette determinant, which describe the geometry of 
geodesic deviation. But this is precisely the property of the geometry encoded in the 
Penrose limit, in which the full spacetime metric in a tubular neighbourhood of the 
chosen null geodesic 7 (with V = x"" = 0) is approximated by 

ds'^ = 2du dV + CabWx" dx^ . (2.3) 

Eq.(2.3) describes a gravitational plane wave, in Rosen coordinates. An alternative 
form, in terms of Brinkmann coordinates, is 

ds^ = 2du dv - hij{u)z'z^ du^ + dz' dz' , (2.4) 

where the profile function hij{u) = R\ju{u), the components of the Riemann tensor 
which occur in the Jacobi equation describing geodesic deviation. An elegant and 
comprehensive review of the Penrose limit can be found in ref. [12, 13], which explains 
the scaling by which the plane wave metric (2.4) arises as the leading term in an 
expansion of the full metric in Fermi null coordinates around 7. The relation with 
vacuum polarization and photon propagation is described in full detail in our earlier 
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work, especially ref.[4], which contains a complete review of the geometry and notation 
used here.^ 

While the relation with geodesic deviation and therefore the motivation for using 
the Penrose limit is most readily seen in terms of Brinkmann coordinates, the field- 
theoretic calculations are best expressed using Rosen coordinates. The relation is given 
in terms of a zweibein E^a{u) for the transverse coordinates in (2.3) and (2.4): 

Cabiu) = E\{u)6,,EMu) , (2.5) 

where E'^a is obtained by solving the differential equation^ 

Eia + hijE\ = . (2.6) 

These equations are solved subject to the requirement that Qij = EiaEj"' is symmetric. 
The coordinates are related by 



z' = E\{u)x'' , 
v = V + kEiJu)E\(u)x''x^ 



(2.7) 



with u common to both sets. The zweibein is not unique and this property is shared 
by the Rosen coordinates which are specially adapted to describing the null congruence 
around 7. 

The Van Vleck-Morette determinant is defined from the geodesic interval 

1 r 



a(x, x'] 



dTg^,{x)x^'x^ , (2.^ 



where x^ = x^{t) is the geodesic joining x = x{0) and x' = x{l) and is given by 

A(x, x') = ^ ^ det . (2.9) 

For a plane wave, with g = g{u), the VVM determinant is a function only of the null 
coordinate u, that is A = A(u, u'). 



^Except that here we are following ref.[5] and using and for the transverse coordinates in 
Rosen and Brinkmann respectively, rather than y° and j/*, and have changed the sign convention for 
V, V. 

^Note that i, j, . . . indices are raised and lowered using Sij while a,b, . . . are raised and lowered with 

Cab- 



-8- 



It is convenient when describing the mode functions on plane waves later in (2.20) 
to introduce the 2x2 matrix %p°'^{u) as the indefinite integral 

^<'\u) = / du Ciu)~^ , (2.10) 



together with ■i/j"'''{u,u') = ip"'''{u) — ip'^^^u'). This is related to the VVM determinant 
by 

A(u,u') = , • 1 ,/ X • (2.11) 

^ ' ' ^Jg{u)g{u') det^(M,M') ^ ^ 

The geodesic interval for a plane wave spacetime then takes the form [3] 

a{x,x') = {u - u'){V - V) + --—^{u,u')-bix - x'Yix - x')' . (2.12) 



An equivalent way to construct and interpret the VVM determinant is by thinking 
of the behaviour of a spray of geodesies which emanate from the point {u', f = 0, = 0) 
in Brinkmann coordinates. In the space-like directions the geodesies satisfy 

^^+h,,{u)z^=0. (2.13) 

If the initial conditions of the spray are chosen as z'^{u') = and dz'^ {u') / du = 6ij then 
we write the solution as z^{u) = Aij{u,u'). The VVM matrix is then simply 

A,jiu,u') = iu-u')[A-\u,u%, , (2.14) 

whose determinant is A{u,u'). 

We will be especially interested later (see section 6) in the special class of sym- 
metric plane waves, or Cahen-Wallach spaces. Here, the profile function hij{u) in the 
Brinkmann metric (2.4) is M-independent. It is convenient to choose transverse coordi- 
nates such that hij = crfSij is diagonal. Then, in Rosen coordinates 

Cabiu) = COS^{(JaU + Ca)6ab ■ (2.15) 

Notice that at can either be real or purely imaginary and Ca are constants.'^ The VVM 



•^If one wants to impose the null energy condition then at least one of the ai is real, say cti, and 
futhermore cti > \(T2\- However, since we are taking the background metric to be fixed and non- 
dynamical there is not really a physical reason to impose this condition. Later we shall look at an 
example with both ctj imaginary since it leads to a particularly simple analytic structure. In fact the 
case with cti = CT2 and both imaginary arises as the Penrose limit of the product geometry of three- 
dimensional de Sitter space with the real line, or circle. This is relevant to a version of the AdS/CFT 
correspondence [29]. 
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determinant for a symmetric plane wave is 



aj(u — u'^ 



AK.O = n-TT^^- (2.16) 

i=l ^ 

Symmetric plane waves are one class of homogeneous plane waves. Another is 
the class of singular homogeneous plane waves, where the profile function is hij = 
^(l — a^)6ij. Like the symmetric plane waves, which have the translation symmetry 
■u — > M + c, the singular homogeneous plane waves also have an enhanced symmetry, 
since the metric is invariant under the scaling u — )■ Am, v — )■ X~^v. As a result, the 
VVM determinant A{u,u') is a function of the single variable u'/u. Explicitly, 

2 



1=1 



This class of plane waves arises as the Penrose limit of certain spacetimes with sin- 
gularities, notably black holes and some Robertson- Walker spacetimes, and will be 
considered in section 7. 



2.2 QFT on a plane wave spacetime 

To study QFT on a background spacetime, we would normally construct an appropriate 
basis of mode functions describing propagation from an initial Cauchy surface. Plane 
waves, however, do not admit Cauchy surfaces. Nevertheless, as described in ref.[26], 
we can still set up an initial value problem in a light-front formalism, choosing the null 
coordinate u to play the role of the time coordinate in the conventional case. While 
the surfaces u = constant are not genuine Cauchy surfaces the initial value problem 
is still well defined when suitable boundary conditions in the transverse directions are 
specified. 

A special feature of plane waves spacetimes is that the solutions of the massive 
Klein-Gordon equation 

(□-m2)0(x) = , (2.18) 
are WKB exact [26]. We can write a basis for the solutions in Rosen coordinates as 

2 

0(x) = e-^$p(x) , (2-19) 
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in terms of the modes 



g{u) ^'''^exp 



(2.20) 



where 'il)°'^{u) is defined in (2.10). The symmetries of the metric allow the definition of 
a 3-momentum p = {u,pa), where a; is a null component and the pa are spacelike. 

The modes are classified as having positive or negative frequency with respect to 
evolution in the coordinate u according to whether a; > or a; < 0, respectively. On 
the null surfaces of constant u, there is a Klein-Gordon inner-product 



($, $') = -i J dV £x ^fg{u) $'*SJ<I> 



with respect to which the modes satisfy the orthonormality property 

{^p,^,)=2uj{27^f5{uj-uj')5^^\pa-p':) . 



(2.21) 



(2.22) 



At this point it is important to recognise, as shown by Gibbons [26], that for plane 
waves there is no particle creation via the usual curved spacetime Bogoliubov transfor- 
mation. The reason is that since dy is a Killing vector, a positive/negative frequency 
mode e^*'^^ will remain a positive/negative frequency mode. As a consequence the 
vacuum state at some inital time remains empty of particles at later times. This plays 
an important role in the next section in simplifying the application of the Schwinger- 
Keldysh analysis. It also means that the type of particle creation arising through a 
non-trivial Bogoliubov relation between in and out vacua cannot be the explanation 
for the increase in amplitude of the photon field observed in some backgrounds. 

The Green functions of (2.18) play an important role in the analysis that follows. 
As usual, we define the Wightman functions 

G+(x,x) = (O|0(x)0(x')|O) , G-{x,x') = (O|0(x)0(x)|O) , (2.23) 

with mode expansions: 

G± = ±^ J ^0{±co) J <fp %{x)%{xye-^2^ . (2.24) 

"^Notice that in general these modes have singularities at points where the metric is degenerate. For 
example, for a Cahen-Wallach space these occur when the cosine functions in (2.15) vanish. However, 
these are only coordinate singularities. 
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The "Feynman" Green function, adapting the definition to our choice of null coordi- 
nates, is then 

iGF{x,x') = (O|r„0(x)0(x')|O) = e{u - u')G+{x,x') + e{u' - u)G4x,x') , (2.25) 

where T„ denotes w-ordering. A similar definition holds for the anti u-ordered Dyson 
function: 

iGY){x,x') = (O|T;0(x)0(x')|O) = e{u' -u)G+{x,x') = e{u-u')G^{x,x') . (2.26) 

Retarded and advanced Green functions, with support in the forward and backward 
lightcones respectively, are also given in terms of the Wightman functions: 

G,et{x, x') = -e{u - u')G{x, x') , Gadv(a;, x') = e{u' - u)G{x, x') , (2.27) 

where the Pauli- Jordan, or Schwinger, function is given by the commutator, 

iG{x,x') = (O|[0(x),0(x)]|O) = G+(x,x') -G_(x,x') . (2.28) 

The Green functions can be written in a "proper-time" representation using the 
expression (2.24) in terms of the mode functions. Performing the Gaussian integral 
over the transverse momentum using (2.10) and (2.12), and using the substitution 
T = ^^^^ in the integration over oo, we find 



GAx,x') = --j^detA{x,x') / ^e^^e— (2.29) 
(47r)2 Jo 12 



for M — m' > 0. Note that, for future use, we have quoted the result for — 2 transverse 
dimensions. Similar results for G+{x,x') for u — u' < 0, and for G-{x,x'), follow from 
the identities G-(x,x') = G+(x',x) = G+(x,x')*. 



The proper-time representations for all the other Green functions follow straight- 
forwardly. For example, 

iGf{x,x') = j.JdeiA{x,x') / ^e^^e-'"^^ , (2.30) 

(477)2 Jo T2 

for all u — u'. 
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3 Light-Front Evolution of Renormalized Quantum Fields 



In this section, we set up and study the initial-value problem describing the evolution 
of a renormalized quantum field in a plane-wave background spacetime, incorporating 
the effect of vacuum polarization. 

3.1 The initial value problem 

The techniques for studying real-time non-equilibrium phenomena in quantum field 
theory were established originally by Schwinger and Keldysh [14-17]. The essential 
idea is to develop an equation of motion, obtained from a suitably defined IPI effective 
action, which describes the time evolution of the "m-m" vacuum expectation value 
j{0in\A\0in) J of the quantum field in the presence of a source, which can be used 
to engineer appropriate initial conditions. In our case, we will consider evolution in 
the light-cone variable u from an initial null surface, but the essential formalism is 
unchanged. 

This should be contrasted with the more familiar IPI effective action derived using 
the conventional Feynman rules, for which the corresponding equation of motion de- 
scribes the 'Hn-ouf VEV j{0out\A\0in) j- The distinction is clearly important when the 
vacua do not coincide, i.e. \Oout) 7^ |0m), and there is a non-trivial Bogoliubov trans- 
formation between the states and associated particle creation. As already emphasised 
[26], this is not the case for plane- wave backgrounds, so we expect our analysis to be 
independent of which formalism we choose. We now verify this explicitly, justifying 
a posteriori the treatment of the refractive index problem in curved spacetime in our 
earlier work. 

The essential features of the Schwinger-Keldysh, or "in-in" , or "closed-time-path" , 
formalism are elegantly summarised in refs. [30-32] and we refer to these papers for 
further details. For simplicity, we describe the formalism first for the QFT of a single 
scalar field 0. The key idea is to let the in vacuum evolve independently under two 
different sources J^(x), comparing in the far future, in the out region. This defines a 
generating functional 



Z[J+,J-] 



expiW[J+,J-] 

J-(0m|0jn)j+ = 



out 



(3.1) 
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The ^Hn-in" expectation value of the quantum field in the presence of a physical 
source J(x) is then given by the derivative with respect to either source, evaluated with 
the sources set equal to J, 



, ,,, V 5W[J+,J- 



5J± 



(3.2) 

j+=j-=j 



The generating functional has a path integral representation in terms of two field 
variables (h^ix) as follows: 



Z[J+,J-] = j T^(t)^ j V<i)- ex^i(s[<p+]- S[<i)-] + J+<p+ - J-(t)-^ 



(3.3) 



where it is understood that the integrals are over field configurations which coincide on 
the out hypersurface corresponding to \ipout)- Evaluating in the free theory gives the 
expression [30-32] 

Z[J+, J-] =exp-^ j d'^x^ j d^x' ^(^J+{x)Gf{x,x')J^{x') 

+ iJ+{x)G.{x,x')J-{x') + iJ-{x)G+{x,x')J+{x') (3-4) 
+ J~ {x)Gj:,{x, x')J^ {x') 



so different derivatives of Z with respect to the two sources yield the full set Gp, Gd, 
G+ and of Green functions. 

Expressions (3.3) and (3.4) can then be used in perturbation theory as the basis 
of an extended diagrammatic expansion, with propagators G++ = —Gp, G+- = iG-, 

G ^ = iG+, G = —Gb linking vertices which are "all +" or "all — ", with opposite 

signs of the coupling. In turn, this allows the construction of a IPI effective action 
r[0+,0~] by a Legendre transform in the usual way, written in terms of fields defined 

U A± _j_SW[J+,J-] 

For our problem, we initially consider the equation of motion for a massless scalar 
field A interacting via an eA0^ interaction with a scalar field of mass m with Green 
functions described in section 2. It is straightforward to see that the required equation 
for the in-in expectation value A{x) = j {Oin\A{x)\Oin) j is 
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(or similarly with and A interchanged), where r[y4+, A ] is the one- loop effective 
action with quadratic part 

T[A+,A-] = -\j d''x^[A+{x)UA+{x) - A-{x)UA-{x)^ 

- d'x^g I rfVV?(A+(x)Gp(x,x')^A+(x') ^3^^ 

+ A+(x)G_ (x, x' fA- (x') + A- {x)G+{x, x'fA+{x') 
+ A-GB{x,xyA-{x')'^ . 

This gives an equation of the form (2.1), with the vacuum polarization 

T^sk{x,x') = '-^{GFix,x'f + G4x,x'f) . (3.7) 

From the definitions (2.23), (2.25) above, we now see that 

UsK{x,x') = --^9{u-u'){G+{x,x'f-G4x,x'f) , (3-8) 

with support only for u' < u. Moreover, given that the commutator [0(x), 0(a;')] is a 
c-number, it follows readily that 

nsK(a:,a:') = -!^^(,, - ^,')(O™|[0(x)^ 0(x')']|O,„) . (3.9) 

Since the commutator vanishes for spacelike separated points, it follows that nsK(a;, x') 
vanishes for x outside the forward light-cone of x' . This confirms that (2.1), with 
the vacuum polarization taken as nsK(a;,x'), is a causal, albeit non-local, equation of 
motion for the in-in VEV A{x) of the field. In contrast, the usual in-out formalism 
with the Feynman vacuum polarization 

nF(x,x') = ^Gf(x,x')' = --^(e{u-u')G+{x,x'f + e{u' -u)G^{x,x'f'^ , 

(3.10) 

is not manifestly causal. 

Now, for propagation in a plane-wave background, the only non-trivial behaviour 
is in the u direction and we can expand the solutions in terms of the basis (2.20) of 
on-shell modes. Accordingly, for a massless field A(x), we take 



A{x) = A{u)^p{x) 



(3.11) 



and choose a; > 0, so that it is a positive frequency solution. So for the simplest case 
where the transverse momentum Pa = 0, we just have (see (3.2)) 



A{x) = yi(M)$(^,o)(x) = A{u)g{ur'/^e'^'' . (3.12) 

Notice that the amplitude factor A{u) here is allowed to be complex and the physical 
solution for A{x) is the real part of (3.12). 

A key result now is that if A{x) has postive frequency, it follows from the form of 
(2.24) that 

dV'G-{x,x'fA{x') = , (3.13) 

for a; > 0. This is a manifestation of cu-conservation at the vertex, in itself a consequence 
of the dy isometry of the plane- wave background. Hence, for solutions of the kind (3.11) 
we can replace IIsk in (2.1) by the Feynman vacuum polarization IIf if we wish, the 
difference vanishing by virtue of (3.13). This is the technical mechanism by which the 
isometry of the plane-wave background, which ensures the absence of particle creation 
and the identity of the in and out vacua, guarantees we recover the same final results 
for "photon" propagation and the refractive index whether we use the conventional 
Feynman or Schwinger-Keldysh formalisms. 

If we now exploit the symmetries of the plane- wave background to define the "par- 
tial Fourier transform" of the vacuum polarization with respect to the modes (2.20), we 
can reduce the equation of motion to an effective one- dimensional problem. Defining 



(3.14) 



(27r)^<5(c^ - io')6^'\pa-p'J IlsK{u,u';uj,Pa) 

= j d^Xy/^ J d^x'^/^ <!?l{x)UsKix,x')<!?p,{x') . 

we can rewrite (2.1) in the form 

{-2iujdu + M'^)A{u) + [ dv! ti{u, u'; oo, p)A{u') = . (3. 15) 



'no 



where for generality we have now introduced a mass M for the A field. For later use, 
we have assumed here that the coupling constant is "turned on" at some specified 
value Uq. This will clarify the later discussion of initial value conditions and the optical 
theorem; we can of course set Mq — ?■ — c>o at any point. The upper limit follows from 
the constraint u — u' > from (3.8). 
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We can further simplify the expression for Il{u,u']u,p) by exploiting translation 
invariance in the transverse space to write it in terms of an integral over the relative 
coordinates, V = V ~ V' and x"" = x"" — x'"'. This gives simply: 




(3.16) 



We now come to the choice of initial conditions. The simplest choice, which we 
call "Type I" , is to assume the coupling is "switched on" at some initial value surface 
u = (so that Uq = 0) and specify the value of the field ^1(0) on that surface. In this 
case, the equation of motion describing the evolution of the field A{u) is simply 

pu 

-2iujA{u) + M'^A{u) + duIlsK{u,u;uj,p)A{u)=0. (3.17) 

Jo 

Prior to switching on the coupling, the field simply evolves according to the tree-level 
equation as A{u) = A{0) exp — ^^^j^. 

With this choice, the field on the initial value surface m = is a bare field, and as 
we follow the evolution for m > it will become dressed with a vacuum polarization 
cloud of virtual (p pairs according to (3.17). For Acj)"^ theory in d = 4, where the only 
UV divergence is absorbed by a mass renormalization, we can follow this real-time field 
renormalization explicitly. We expect an initial fall in A{u) as the bare field becomes 
screened by the dressing. 

As with transient phenomena in general in quantum field theory, this instantaneous 
"switching-on" of the coupling may be considered rather unphysical. However, we can 
simulate these initial conditions while allowing the coupling to be non-vanishing for 
all u {uq — )■ —oo) by introducing an appropriate source J{u) on the right-hand side 
of (2.1) chosen to cancel the contribution of the integral for — oo < m' < 0, recovering 
(3.17). However, this requires J{u) to be non- vanishing also for u > and it is arguable 
whether the introduction of such a fine-tuned source for all u is really any more physical 
than the original model of switching on the coupling. As we will discover, studying 
the initial transient dressing even in flat space provides an important insight into the 
long-time behaviour of the field as it evolves in time-dependent curved spacetimes. 

In the second part of this paper, we study examples of field theories with different 
short-distance behaviour, specifically A(f)'^ in d = 6, which requires a UV divergent field 
renormalization but which is asymptotically free, and QED, which is IR free and where 
perturbation theory in the short-time, UV regime is afflicted by the Landau pole. In 
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these cases, we find a variety of problems witli tlie Type I initial conditions. Since 
our primary aim is to study tlie universal long-time behaviour of the field, which we 
would expect to be independent of the preparation of the initial state, we also consider 
alternative initial conditions. 

In "Type II" initial conditions, which have previously been used to study relaxation 
problems in non-equilibrium QFT (see, e.g. ref.[18]), we use the source to hold the field 
A{u) fixed from mq — )■ — oo to the initial surface u = 0. In this case, the equation of 
motion is: 

{u, u'; uj,p)A{u') 

Jo 

+ 1 du'tl^K{u,u-u,p)A{Q) = Q ^^-^^^ 

J —oo 



In this scenario, the field at the initial value surface m = is already renormalized and 
partially dressed, and as m > we watch it relax from this state. This choice of initial 
conditions circumvents the problems with short-distance physics for theories such as 
QED, where the Type I conditions are not controllable in perturbation theory. 

A further variant, "Type III", is to constrain the field to evolve for m < with a 
phase yi(0) exp — where is the renormalized mass (to be distinguished from 
the bare mass in the second term in (3.18)). The new equation of motion is the obvious 
generalisation of (3.18). In this case, the field at m = is already renormalized and fully 
dressed and, as we confirm in section 6, no further evolution occurs in fiat spacetime. 
However, in curved spacetime, these initial conditions are particularly appropriate and 
allow us to study the effects of curvature on a fully-dressed quantum field, even in cases 
requiring a divergent UV field renormalization. 



3.2 Vacuum polarization 



We now sketch the evaluation of the vacuum polarization in a plane-wave background 
in the form we need here. Further details, including the equivalent calculation for QED, 
can be found in refs.[3, 4]. We consider initially the Feynman form, 

Ii^{x,x') = --^G^{x,x'f . (3.19) 

First, consider briefly the equivalent calculation in flat spacetime. Here, a{x,x') = 
|(x — x')^ and A{x,x') = 1, and because of translational invariance the self-energy only 
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depends on the relative position. We can therefore Fourier transform with respect to 



the four relative coordinates x'^ 



x'^. Writing the two proper-time variables as 



Ti and T2, we then change variables from (Ti,T2) to (T, ,^), where T = Ti + T2 and 



dTi dT2 



°° dT di 



and find: 



d xe 



ip-x 



U(x) 



2(47r) 



d^ 



dT 



-i(m2+p25(l-f))T 



(3.20) 



(3.21) 



This displays a UV logarithmic divergence as the proper time T — )■ 0. This is a 
conventional UV divergence that is cancelled by a mass renormalization for the A field. 
We can evaluate the integral by Wick rotating the proper-time T — t- —iT and by 
introducing an explicit cut-off A^^ on the T integral, where A is a momentum scale: 



2 (47r)- J, 



di 



00 jrp 

T 



A-2 



2 (47r) 



d^ log 



m 



A2e-7i5 



(3.22) 



This is the standard result for the one-loop contribution in momentum space. 

Now consider the calculation of the vacuum polarization in curved spacetime from 
(3.16), initially using the Feynman form Ilp{x,x'). The Green functions are given in 
(2.30). The integrals over f are Gaussian, while the integral over V generates a delta 
function: 



dV exp 
d'^x exp 



4rg(i-g) - u')ilj{u,u')-b^''£'' exp[-ipax°] 



(3.23) 



4?7rT^(l -O-^St 77^ 



[u — u 



[u — u) 



Since u — u' = 2c<;^(l — ^)T, the conditions u > and T > mean that u > u'. This 
is the calculational mechanism already encountered in (3.13) which ensures that, when 
acting on a positive frequency solution, the Feynman self-energy is only non- vanishing 
when u > u' and therefore gives the same result as the retarded, Schwinger-Keldysh, 
self-energy. 
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Imposing the delta function constraint, cancelling the p"-dependent terms in (3.23) 
against the identical terms in the mode functions, and collecting terms, we find: 

np(^, vl- 00, p) = -^j^,e{u - u') ^ > / d^e-^-^i^-O . (3.24) 

2 (47r)^ u — u' Jq 

The usual UV divergence now appears as the singularity in the limit u ^ u' and 
can be renormalised in the standard way. Crucially, there are further "geometric" 
singularities when u and u' are conjugate points on the null geodesic describing the 
classical "photon" trajectory. At these points, the VVM determinant diverges. As 
explored in detail in refs.[3], it is these singularities which give rise to the novel analytic 
structure of the refractive index in curved spacetime which modifies the conventional 
form of the Kramers- Kronig dispersion relation while maintaining consistency with 
causality. The correct prescription for dealing with these singularities follows from the 
Feynman prescription in real space; namely u — u' u — u' — iO^. 

Alternatively, we can derive the Schwinger-Keldysh vacuum polarization function 
IIsk{u,u'',<^,p) directly from (3.16) with IIsk{x,x') — )■ —ie^6{u — u')G^{x,x'Y, imme- 
diately exploiting (3.13). Writing the Green functions directly in terms of the modes 
according to (2.24), we have: 

Usk{u,u';u,p)= --9iu-u'W9{u)9{u') j d'^ J J ^^^^^^ 

X e--^("-"')(^ +^)$(,)(x)t$,,(x)<l>,,(x')^$P2(x)$p2(a:')^<f(p)(^') , 

where d^pi = dui dp} dpj and the integrals over uji are restricted to cjj > 0. The integrals 
over the relative transverse positions x impose momentum conservation 

v"=Vi+V2, uj = uji+uj2. (3.26) 

It is convenient to solve the second constraint by taking ui = uC, and U2 = {1 — 
with < ^ < 1. This leaves a Gaussian integral over p^: 

^e^^, ,x r / N / /Ni-1/4 ,^ 1 ""^("-"') 



ie -1/4 I 1 

nsK(M,M';w,p) = - —B[u - u!)\g[u)g[u:)\ J ^^^^^ _ ^^ e" 

d'^Pi 



exp 



(2vr)3 n 2uai-0 



(Pl - ipTi^ab{u, U){pi - ^p) 



= -oTT^^ / / die , 

2 (47r)^ u — u' Jq 

(3.27) 
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where we have again used (2.11) to write the result in terms of the VVM determinant. 
Naturally, this reproduces precisely the expression (3.24). 



3.3 Field evolution and the dynamical renormalization group 

The solution to the initial value problem, for Type 1 initial conditions, is given by 
solving the equation of motion (3.17) for the u-dependent amplitude factor A{u) defined 
in (3.11). We write the solution as a perturbative expansion in as:^ 

A{u) = Aiuo)e- 2^ (1 + iQ^^\u) + • • ■ ) . (3-28) 
Implementing the constraint that IIsk{x,x') has support only for u > u', we find: 

Q^^\u) = -7r / du'IisK{u",u'-u,p)e 2;:^ , (3.29) 



SO from (3.27), 

Re Q^^\u) is of course UV divergent and, for the c? = 4 Acj)"^ theory, this divergence is 
absorbed as usual in a mass renormalization of M^, as shown explicitly in section 5. 

The imaginary part, lmQ^^\u), controls the amplification or attenuation of the 
amplitude and plays a vital role in our analysis from now on. Using the symmetry 
property A{x,x') = A{x',x), we can show that this is given by an expression almost 
identical to (3.30), but with the integration range extended: 



(3.31) 



This expression, with the integration limits shown, is key to clarifying issues related 
to positivity and the optical theorem in the next section. Note also that Im Q^^\u) is 



^For Type II or III initial conditions, with initial value surface u = 0, the solution for w > is 
written as 

A{u) = yi(0)e"^i:r (1 + iQ^^){u) - iQ(i)(o) + •••). 

Note that at this order, the equations of motion for Type III and Type I initial conditions are identical 
for u > 0. 
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free of UV divergences since the u' integral avoids the singularity at u' = u" by virtue 
of the Feynman prescription u" — «'—)■ u" — u' — iO"*". 

Now, for backgrounds which are translation invariant in u (the symmetric, or 
Cahen-Wallach, plane waves), the VVM determinant will be a function only of the 
separation of the two points, i.e. A = A(x — a;'). It follows that for large u, Q^^\u) ~ u. 
When u is large, therefore, the 0{e^) correction will itself be large and perturbation 
theory should break down. In fact, such large secular terms are a common occurrence 
in perturbation theory and in this translation-invariant case, the problem is overcome 
by a resummation which exponentiates the dependence on Q!^^\u) in (3.28). 

More generally, this kind of resummation is performed by what is known as the 
dynamical renormalization group (see [18] and references therein). This is a way of 
absorbing the secular terms into a renormalization of the amplitude A{u) at an arbitrary 
time u*. The quickest way to arrive at the re-summed formula is to write in perturbation 
theory 

A{u) = A{u*)e-'^^^^{1 + z[Q«(«) - Q^'\u*)] + ■■■). (3.32) 
Since the right-hand side cannot depend on the arbitrary scale u* we have, to O(e^), 

= —A{u*) + i A{u*)-i , ^ ' A{u*) . (3.33) 

du* 2uj du* 

Solving this, using the initial condition Q^^\uq) = 0, and finally setting u* ^ u gives 

A{u) = A{uo)e' 2:^ e*2(i)H _ (3.34) 

The effect of the resummation is to exponentiate the first order perturbative correction. 
The importance of the dynamical renormalization group is that we can apply it even 
in the non-translationally invariant cases. 

Notice now that a positive imaginary part of Q^^\u) signals an exponential decay 
of the amplitude, which is the standard expectation. For translation-invariant back- 
grounds, where lmQ^^\u) is proportonal to u in the long-time limit, this allows us 
to identify a constant decay rate. However, this interpretation does not extend to a 
general non-symmetric plane-wave background where the usual concept of a decay rate 
is not well-defined. 
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3.4 Refractive index and analyticity 



At this point, we make contact witli our previous work [3, 4] on photon propagation 
and the refractive index for curved spacetime. 

In the massless (M = 0) case, after applying the dyanmical renormahzation group 
resummation, we have found the solution for the field A{x) in the form 

A{x) = ^K)$p(x)e'2*'^(") , (3.35) 

where 

= du'llsK{u,u : 00, p) . (3.36) 



Note again that whereas previously we used the Feynman, in- out, formalism, all our pre- 
vious results for the refractive index are identical with those found using the Schwinger- 
Keldysh, in-in, formalism. Comparing with refs.[3, 4], we identify the refractive index 
associated with the solution (3.35) as 

n{u;uj) = l + -^Q^'\u) (3.37) 



and so 



n{u]uj) = 1- —J I du'IlsK{u,u';uj,p) . (3.38) 

Juo 



To complete the identification of notation, substitute the explicit expression (3.30) 
for Q^^^{u) and introduce the variable t = u — u' . We then have 

1 r""" dt ^m^t , 

n{u- u) = l + di - e ^"^(-o ^/\{u,u-t) . (3.39) 

The upper limit of the t integral goes to oo, as in our previous discussions of the 
refractive index, when we take the switch-on time Mq — )■ — oo. Of course, this expression 
still requires renormahzation. This will be considered in this formalism in section 7, 
after the detailed analysis of translation invariant spacetimes in sections 5 and 6 using 
Laplace transform methods. In ref.[3], we also introduced a variable z = 2u]l(i-^) 
wrote (3.39) in the form: 

n{u; cu) = l + j^-^ [ d^ ^(u; z) , (3.40) 
(47r)Mcu^ Jq 
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with 



J^{u;z)= -e-''' ^/A{u,u-t) . (3.41) 
Jo t 

Apart from the obvious differences because here we are considering a pure scalar theory 
rather than QED, these expressions match the forms presented in [3, 4], in particular 
section 6 of [3]. 

From this point, ref.[3] discusses in detail how the singularities in the t-plane, which 
arise from the VVM determinant when u and u — t are conjugate points along the null 
geodesic describing the classical photon trajectory, integrate up to give a novel analytic 
structure in the complex w-plane for the refractive index. In particular, the presence 
of new cuts related to the background geometry violate hermitian analyticity, which 
is a standard property of S-matrix theory. The resulting modification of standard 
theorems of flat-space quantum field theory, especially the Kramers- Kronig dispersion 
relation, allow the apparent paradoxes associated with low-frequency super-luminal 
phase velocities - that is, R.en{u; 0) < 1 - to be resolved while causality is maintained. 
In what follows, our emphasis is different as we try to reconcile the amplification of the 
field implied by a negative imaginary part of the refractive index with unitarity. 



3.5 Dressing and undressing a quantum field 



One of the most remarkable features to emerge from our analysis of the refractive index 
for QED was the discovery of cases where Imn(M;a;) is negative, corresponding to an 
amplification of the field as it propagates through the background curved spacetime. 
To begin to develop some intuition into this phenomenon, we consider first the weak 
curvature limit. Here, we can expand the VVM determinant to linear order in the 
curvature 

A{u,u') = 1 + ]:Ruu{u){u - u'f - 1-R^^{u){u - u'f + ■■■ (3.42) 
12 

where the Ricci tensor R^u = R\iu- Substituting this expansion into the integral (3.36), 
and after UV regularization, we determine the correction to the refractive index: 

-) = 1 - - (i^^^"-(^) + ■ ■ • • ^3.43) 

The analogues of this expansion for scalar and spinor QED appear in appendix A and 
(1.1) for the latter. In this approximation, the ratio of the amplitude for two points 
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Ml, U2 on the null trajectory 7, in the long-time limit, is 



A{ui) = A{u2) exp 



■{Ruu{Ui) - Ruu{u2)) 



(3.44) 



_(47r)2 840m6 

If the Ricci curvature component Ruuiu) decreases along 7, then from the Jacobi equa- 
tion (see [4]) we see that the rate of acceleration of the nearby geodesies is positive, 
i.e. the tidal forces are increasing in the direction of stretching the virtual cloud. In 
this case, the amplitude decreases because more virtual pairs are being produced: A 
is becoming more dressed. On the contrary, if Ruu{u) increases along 7, i.e. the tidal 
forces are increasing in the direction of squeezing, then the amplitude increases because 
virtual pairs are recombining: A is being "undressed" and is returning to its bare state 
[5].^ These effects are illustrated heuristically in Fig. 1 with the role of the photon 
being played by A and the e^e~ pairs by cj) pairs. 

We now come to the second, non-perturbative, mechanism which can produce an 
imaginary part for the refractive index. If Ruu{u) is constant, in which case the metric 
(2.4) describes a symmetric plane wave, or Cahen-Wallach space [19], then to linear 
order in the curvature the amplitude is constant. In fact, lmn{Lj) = to all orders in 
the curvature expansion. However, in the case where at least one of the eigenvalues of 
the constant hij is negative, there is a positive contribution to Im?T,(a;), so a decaying 
amplitude, which is non-perturbative in the curvature. The simplest case to consider 
is hij = —a'^6ij, which is discussed in detail in section 6. In this case, the VVM 
determinant is 

2 



u 



we then have 



1 



sinh aiu 



dt 



u' 



niu] 



die 



(3.45) 



(3.46) 



{AtxYAuj'^Jq sinh (at) jo 
The imaginary part is then given by (3.46) by extending the contour to run from —00 
to +00 just under the real axis. \m.n{uj) can then be computed by deforming the 
integration contour so that it picks up the residues of all the poles along the negative 
imaginary axis at t = n = 1, 2, . . .. This gives 

^2 _ nl 



TT 



(47r)2w2 



di 







(3.47) 



related physical picture for particle creation in a background gravitational field, based on quan- 
tum mechanical tunnelling, has recently been presented in ref.[33]. This shares our basic intuition of 
relating these curved-spacetime phenomena to geodesic deviation through the Jacobi equation, though 
here we are considering one-loop quantum field theoretic effects due to vacuum polarization. 
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which is non-perturbative in the curvature. In this case the result is independent of 
u and so it imphes a constant rate of production of pairs. This process would be 
kinematically disallowed in flat space, but there is no such threshold constraint in curved 
spacetime. This kind of below-threshold decay is a well-known property of curved space 
arising from the lack of time-translation invariance. In particular the situation in de 
Sitter space has been known about for a long time [36-39]. 

This result essentially confirms the proposal of ref. [28] that interacting quantum 
fields in a plane wave background with negative eigenvalues of hij will decay. We have 
shown that the process is perfectly consistent with unitarity and if one starts from an 
initial value surface then only a finite number of 0^ pairs will be produced per unit 
volume. In the next section we turn to the issue of fate of the optical theorem in curved 
space. 

4 Optical Theorem in Curved Spacetime 

The observation that the imaginary part of the refractive index can be negative, at 
least for trajectories and spacetimes where the metric in the Penrose plane- wave limit 
is M-dependent, brings us to the critical issue - how can we reconcile the increase in 
the amplitude associated with lmn{u; u) < with unitarity and the optical theorem? 

The optical theorem is familiar from the discussion of scattering amplitude in flat 
space: in the context of the present model, it states that the imaginary part of the 
one-loop vacuum polarisation diagram yields the rate for the tree-level process A — )• 00 
and is therefore non-vanishing only above threshold. In curved spacetime we need a 
reformulation of the theorem that avoids the use of asymptotic paticle states. The 
initial-value problem provides a suitable formalism. For clarity, we consider Type I 
initial conditions where the interaction is switched on at m = mq. We then calculate 
the total probability for the tree-level process A — )■ 00 to occur between times uq and 
u and compare with the solution of the initial value problem for A[x). This gives a 
generalised version of the optical theorem that holds in curved spacetime and directly 
relates the dissipation of the A field to the creation of pairs. 

The resolution of the apparent paradox that lmn{u] u) can be negative lies in the 
fact that we have only shown that the amplitude can increase over a local region. In 
order to verify unitarity, and relate Im 11 to the (positive) probability of production of 
pairs, we have to integrate over the whole region from the "switch-on" surface uq to 
u (see eq.(4.2) below). It turns out that while the amphtude A{u) can increase locally. 
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compared with its value at the initial value surface it must always decrease. The initial 
transient dressing of the A field plays an important role since it is only a field that is 
already dressed that can increase its amplitude by "undressing" ; a bare field, however, 
cannot be undressed any further. 

To derive the optical theorem, note first that 

/ du"\mn{u"]uo) = -\mQ^^\u) = du" du' Imtls^iu" ,u']UJ,p) . 

Juo ^ Juo J no 

(4.1) 

where \mQ^^\u) is given in eq.(3.31). We also leave M = 0, so it is appropriate to 
use the refractive index nomenclature as in section 3.4. Otherwise, IisK{u, u']UJ^p) in 
(4. 6), (4. 8) is accompanied by a simple M-dependent phase, as in (3.34). 



Now consider the the transition probability for the tree-level process A — )■ 
[ d^pi d^p2 



M J (27r)32a;i {27t)^2uj2 



d'^x a/^ A^^^ (x) {x)4'p2 i^] 



UO 



(4.2) 



Here, (f)p{x) are the massive on-shell modes defined in (2.19) and the final-state phase 
space integrals are over the three-dimensional p = {uj,Pa). The normalization factor 
J\f = {A^^\ A^^^) = 2co5^^\0) cancels against an overall integration on {V,x"'). All the 
frequencies, u, ui and U2 in (4.2) are positive since they refer to physical particles. 

The evaluation of (4.2) follows closely the calculation of IIsk{u,u';u,p) directly 
from the mode functions in eqs.(3.25) -(3.27). For convenience, recall here the definition 
of nsK(a;,a;'): 

nsK{x,x') = --^9{u~u')(G+{x,x'f-G4x,x'f) , (4.3) 
the expression for the Wightman functions in terms of the modes (f)p{x): 

G±{x, x') = ±l 9{±co) M^)M^'y ■ (4-4) 

and the observation (3.13) that / dV'G_{x, x'^A^^^^x') = for positive frequency u. 
It then follows that 



e2 



MO 



Pa-,^^{u) = jj I du"^ I du'^ I dV / d'x'<^p{x"yG+{x",x'f%{x) . 



(4.5) 
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Separating out an overall volume factor using translation invariance in the transverse 
space, as in (3.16), and using the properties G+(a;,a;')* = G-{x,x') and G+{x\x) = 
G_{x, x'), we find: 

Pa^,I>4>{u) = / du" / du lm.ti^K{u" ,u']UJ,p) , (4.6) 

where note again that we could equally have written IIf above, and so 

Pa^4,(I,{u) = Au / du'lmn{u"]u) . (4.7) 

This is the statement of the optical theorem in curved spacetime. Unitarity is respected 
and the integral over the whole trajectory of \m.n{u\uj) is positive. This shows how 
unitarity prevents a local amplification of the amplitude from becoming unbounded.^ 

We can also write a local version, defining T{u) = duPA~^<t><t>{u) as the "instantaneous 
rate" of pair production. Then, 

T{u) = 4:u\m.n{u]uj) = / du' ImnsK(M, u'; w,p) . (4.8) 

^ Juo 

Unlike the integrated form, there is no positivity constraint on (4.8) in general. This 
explains why there is no confiict with unitarity in the examples where we have found 
lm.n{u\uj) < 0. However, in cases where we have translation invariance along the 
photon trajectory, as in fiat spacetime, n{uj) becomes w-independent in the large u 
limit beyond the transient region (equivalently Uq — )■ — oo; see (3.39)) and P4_j.00(m) is 
then proportional to m, so unitarity implies Imn(a;) > and F > can be properly 
interpreted as the rate oi A cfxj). 

To summarise, we have shown how the effect of gravitational tidal forces on vacuum 
polarization can alter the dressing of a photon as it propagates through space. In 



^ To emphasise the importance of integrating over the entire region from the "switch-on' surface 
Mo in these expressions, consider the change AP(u2|i*i) in the total probabihty of pair production 
between two nuU surfaces ui and U2- From (4.6) this is: 

AP(u2|ui) = / du" du' lmilsK{u" ,u']uj,p) . 

^ Jui Jua 

Contrast this with the integral I{u2\ui) defined as AP{u2\ui) but with the lower limit of the u' integral 
set to Ml. This is a positive quantity, since by the above construction it can be written as a | • • • P as 
in (4.2), but it does not represent the probability of pair production between ui and U2. In fact, while 
P{ui), P{u2) and I{u2\ui) are all positive quantities, because of the extra region of integration from 
the switch-on surface uq to mi, lS.P(u2\ui) itself need not be positive. 
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particular, we have seen how the imaginary part lmn{u]u) of the position-dependent 
refractive index can be negative as well as positive, corresponding to "undressing' of 
the photon rather than the conventional dressing. Two mechanisms were identified. 
The first, of order udu^/m'^, admits an intuitive interpretation in terms of curvature 
variations along the A field ("photon") trajectory altering the balance of the bare field 
with its virtual e~^e~ cloud, with increasing stretching (squeezing) giving rise to more 
dressing (undressing). The second, which is non-perturbative in cj^iH/m^, can occur 
even when the curvature is constant and is related to the existence of conjugate points 
on the photon's null geodesic. 

Nevertheless, unitarity is still respected and the optical theorem still holds in curved 
spacetime. In its integrated form (4.6), (4.7), it relates the total probability for e~^e~ 
pair production to the integral of lmn{u;u) along the photon trajectory, which is 
manifestly positive. Except in the special case of translation invariance along the null 
geodesic, the corresponding local form (4.8) has no positivity constraint and the usual 
interpretation of lmn{u] u) as the rate of pair production can break down. However, it 
does describe the variation of the amplitude and its interpretation in terms of dressing 
and undressing of the photon field. In a sense, photon propagation in curved spacetime 
resembles the initial transient phase in fiat spacetime, with the characteristic features 
of an oscillating amplitude and below-threshold decay. 

5 Real-time Field Renormalization 

In the second part of this paper, we illustrate these general principles with a range of 
examples featuring different quantum field theories, initial conditions and plane-wave 
backgrounds. 

Most of the discussion will focus on spacetimes with w-translation invariance. In 
this particularly powerful method of analysing the equation of motion for A{u) 

is provided by the Laplace transform, which automatically re-sums chains of one-loop 
vacuum polarization diagrams and provides a convenient formalism in which to exploit 
their analytic structure. 

After introducing the Laplace transform method, the remainder of this section 
specialises to flat spacetime. As well as providing a relatively simple example to develop 
our techniques, it turns out that by studying real-time field renormalization in the 
transient region in fiat spacetime, we develop a lot of insight into the evolution of a 
renormalized quantum field in a curved spacetime background. 
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5.1 Laplace Transform and w-translation invariance 



To introduce the Laplace transform method [18], consider first the ^40^ theory in d = 4. 
To simplify notation for these u-translation invariant cases, we define the kernel function 
in the equation of motion as = Il{u,u'; uj,p) and note that = 

where as in section 3.4 we let t = m — u'. 

The equation of motion, with Type 1 boundary conditions and = 0, is then (see 
(3.17)): 

-2iuA{u) + M^A{u) + dtJ:{t)A{u-t) = . (5.1) 

Jo 

Since the final term is a convolution, this can be solved by introducing the Laplace 
transform, 

/•oo /"OO 

A{s)= due"''A{u) , t{s)= dte-''i:{t) . (5.2) 
Jo Jo 

The transform of (5.1) is simply 

-2iuj{sA{s) - A{0)) + M^A{s) + t{s)A{s) = (5.3) 

and so 

A{s) = . 5.4 

^ ^ 2tus - M2 - S(s) ^ ' 

The inverse transform is then implemented by the usual Bromwich integral 

A(u)^n^ . (5,5) 

J,_,^ 2m 2iujs - M2 - S(s) 

where the contour lies to the right of any singularities in the complex s plane. 

Now consider the renormalization of the theory. The analytic structure for the 
inverse transform depends on properties of the vacuum polarization kernel S(s) and 
is mo del- dependent. Suppose first (this will be the below threshold case in fiat space) 
that there is an isolated pole at s = Sq given from (5.4) as the solution of 

2iujSQ - M| - Sij(so) = , (5.6) 

where we have explicitly indicated the bare quantities in (5.5) with a subscript. We 
renormalize by identifying the mass counterterm, M\ = + SM"^, in such a way that 
M is the physical mass defined as the position of the single particle pole, i.e. 

2iujSo = . (5.7) 
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It follows from (5.6) that the counterterm is 

.M=^-E«(|!). (5.8) 

which is real if we are below threshold. The expression A{u) in (5.5) then takes the 

same form with the mass interpreted as the renormalized mass and the kernel 
replaced by the subtracted form 

SW = E,W-E,(|!). (5.9) 
Field renormalization will be considered in section 5.3. 



5.2 in d = 4: transient evolution in fiat spacetime 

In a M-translation invariant spacetime, the Fourier transform of the kernel Xl(t) is simply 
the usual vacuum polarization with momentum p = {u,—is — \CabP°'P^ ,p"')i which 
implies = —2iujs. So in flat spacetime, we can immediately read off the expression 
for the bare S(s) from (3.22): 

^.(.).n.(-«..)4^/'.,io.( '"'-^y4t-a- ) . (5.10) 

Implementing the mass renormalization subtraction (5.9), we therefore have 

The analytic structure of the Laplace transform is illustrated in the left part of 
Fig. 2. S(s) itself has a branch cut at s = — ^Hi!.^ equivalent to = —Aim?, which is 
the threshold for the production of pairs of particles. It is convenient in the following 
discussion using the spectral density to introduce the variable v"^ = —p^ = 2iojs. Notice 
that with this definition, the u plane is a double cover of the s plane. Consequently, 
whereas in the s-plane there is only a single 2-particle cut beginning at the branch-point 
s = —^^5^, in u plane there are two threshold branch points at u = ±2m and we take 
the cuts to lie along — oo < jy < —2m and 2m < u < oo. 

The spectral density, which plays a key role in our analysis, is then given as usual 
by the imaginary part of the vacuum polarization kernel, 

p{u) = --lmt(^] , (5.12) 



TT \2iuj 
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threshold 



2mf 
iuj 



Mi 
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above 
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y • « 


> 



Figure 2. The analytic structure of the Laplace transform in the below thresold (left) and 
above threshold (right) situations. The diagrams show the particle pole at s = ||j and the 
2-particle threshold branch point and associated cut at s = In the case above threshold 
we have deformed the cut to expose the simple pole on the un-physical sheet. 



defined implicitly with u ^ u + iO~^, i.e. as the limit from just above the cut. For 
theory in d = 4, it is related to S by the once-subtracted dispersion relation, 

■2 ^ roo 



\2iuj 



n^)-^ -77-, 



v' v' — V — zO" 



(5.13) 



From the one-loop expression (5.11), we find the imaginary part: 

—Q{v-2m)Ti / di 
^ h- 



ImE 



'2 (47r 



^^(^-2m)A/l 



2 ' 



(5.14) 



where ^± are the roots of — ^ + ^ = 0. The spectral function for theory in 
(i = 4 is therefore: 



(5.15) 



2 ■ ■ V V- 

valid for z/ > 0. The v dependence of this expression explains the need for the subtrac- 
tion in (5.13) to ensure convergence of the integral for large v. For negative we have 
p{—v) = —pijj). Notice that the imaginary part of S(s) and the spectral density are 
UV finite quantities. 



Below threshold M < 2m, the contribution from the pole and cut are disentangled 
and the result is simply a sum of the two contributions. 



A{u) = Apoic{u) + A 



[u 



(5.16) 
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The single particle pole contributes 



■^pole('w) 



A{0)e- 



1 - ^S'(so) 1 



1 X / Am-' 



(5.17) 



(47r)2 2A/2 ^ V M2 

where, using (5.11) and evaluating the ^ integral, we have 

z 1 



: arctan 



(5.18) 



with z 



4m^ 
M2 



So the contribution from the pole is simply the classical solution with a 



quantum corrected amplitude. The contribution from the cut is: 

up{u)e 



2u 



A„Ju) = -2A{0) I du ■ ~ , 2..2 

2m [u^-M^-ReJ:{^)] +7c^p{uy 



(5.19) 



We can immediately check consistency by evaluating the sum rule^ 



1 + 



1 ^/Am^ 



(47r)2 2M2'' VM2 



2 / dv . , ., 



1 . 



(5.20) 



It is also important to notice that the solution that we have found using the Laplace 
transform remains perturbative throughout the evolution. In other words, in the inverse 
transform (5.5) we can consistently expand in powers of of the coupling: 



A{u) = A{0) 



c+ioo 



ds 



Ml 

2iu} 



+ 



\^ 2iuj I 



+ ■ 



(5.21) 



The pole at s = yields the terms 

A{0)e~^ 



1 + S'( ) +■ 



(5.22) 



which are precisely the terms of order O(e^) in the expansion of (5.17). At this order, 
the cut contribution is 



^cut(M) =^(0) 



(47r)2 



r] ^ 

2„ \v^-M^Y 



1 ^ e 2^ 



(5.23) 



^This follows from the fact that at w = we can pull the contribution from the pole and the cut 
off to a large circle at infinity. In this limit, 2za;s — — E(s) — >■ liujs and the so the contribution is 
a simple pole at infinity with residue 71(0). 
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below threshold 




>u 

Uo 



Figure 3. The amphtude of the field as a function of u for the below and above threshold cases 
for a representative choice of parameters. The case above threshold illustrates exponential 
decay while the case below threshold oscillates transiently before going asymptotically to a 
constant which gives the finite wave-function renormalization . 



Evaluating this integral and comparing with /(^p 
(5.20) explicitly at 0{e'^). 



we readily verify the sum rule 



Since the spectral density is positive, for positive z/, the contribution to A{u) from 
the cut is positive. Although it cannot be evaluated analytically, we can calculate its 
large u behaviour, which is dominated by the form of the integrand near the branch 
point V = m. This gives the long-time behaviour. 



J\.r. 



up to M-independent prefactors. 



2u 



(5.24) 



The combined result for A{u) is plotted numerically in Fig. 3. For large m, the cut 
contribution to the amplitude itself (factoring out the overall e~*^^ phase), oscillates 
with a scale uq = ^,^2'^^^ while decaying as a power law on the scale ud = 
and ultimately A{u) converges to the constant Apoic{u)- Note that the scale of the 
oscillations increases as approaches (from below) the threshold 4m^. We recognise 
the ratio 



A{oo) 



A{0) 



1 



dll{p) 



Qp2 



-M2 



-1 



1 



1 + 



(47r)2 2M2 J V m2 



(5.25) 



'in the large u limit, we use the fact that du ue 2aj — 2m is given approximately by 



exp I 



2im u 



2 ) 
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as the wave-function renormalization factor evaluated in the equihbrium theory. In 
particular, due to the positivity of the spectral density the contribution from the cut 
is positive and therefore the sum rule (5.20) implies Z < 1. So intuitively what is 
happening is that we are seeing the dressing of the field A in real time by the creation 
of (p pairs even though we are below the threshold for decay. The fact that such a 
decay cannot happen energetically is by-passed here because we are looking over a 
finite region of time and consequently the energy has an associated uncertainty which 
allows the below-threshold process to occur. It is important that this effect does not 
involve a constant rate but is just a transient effect. 

We now discuss what happens above threshold, when M > 2m. In this case the 
position of the single particle pole (5.6) sq moves off the imaginary axis. To C(e^) 

^°=2^ + 2^^-K2^)+---- ^'-''^ 

Notice that the mass counterterm is now specified more precisely as 

SM' = -Ret(—] , (5.27) 



Since the 2-particle cut arises form a square-root branch point, the function S(s) is 
defined on a 2-sheeted cover of the s-plane with the Bromwich contour on what we call 
the upper sheet. A closer analysis reveals that with the 2-particle cut lying along the 
negative imaginary axis in the s-plane, the single particle pole lies on the lower sheet. 
Hence, the Bromwich integral only receives a contribution from the cut. However, in 
order to determine the large u behaviour we can deform the cut to the left of the pole 
as illustrated in the right-hand side of Fig. 2 and then the pole does contribute 



A 



pole 



A(0)e' 



M'^U _ IT 

2iuj 2ui 



p(M)u 



2iu \2iu 



+ ... 



(5.28) 



to order 0{e^). So above threshold, for large u the A decays into </> pairs with a 
characteristic life-time 



— p{M) = - — - — 
2u ' {A-nY^u 



1 - 



Am? 



(5.29) 



The behaviour of the field is illustrated in Fig. 3. Again note that as approaches 
threshold (from above), the lifetime becomes large. 
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5.3 A<p'^ in d = 6: field renormalization and initial conditions 



We now consider the initial value problem for a theory, in c? = 6/*^ which requires 
a UV divergent field (wave-function) renormalization. The theory is asymptotically 
free and short distance physics remains perturbative. Nevertheless, we would expect 
difficulties with applying the Type I initial conditions due to the UV divergence when 
the interaction is switched on instantaneously at the initial value surface u = uq = 0. 
Indeed, this is what happens, and we illustrate these difficulties before moving on to a 
solution of the initial value problem for Type II and Type III initial conditions, where 

Mo — )■ — OO. 

Incorporating a field renormalization to 0{e^) in the equation of motion (3.17), we 
have 

-2iuZA{u) + MlZA{u) + / du'tiB {u,u';u,p)A{u') = . (5.30) 

Jo 

The solution for A{s) is then 

A(s) = ^^-3 . (5.31) 

{2ioJS - Ml)Z - Sb(s) ^ ' 

The Laplace transform kernel Sb(s) is again obtained from the flat space liBip^) 
by the substitution — —2ius. In dimensional regularization, this gives 



- 2ius^{l - 



2p-7 



(5.32) 



for the bare kernel. Physical mass renormalization is as before, 

K{s) = ts{s)-ts[^) , (5.33) 

where M is the renormalized mass, and the field renormalization corresponds to the 
further subtraction 

t{s) = tr{s) -{Z-1) {2iujs - M^) , (5.34) 



""^"The extra dimension for d ~ 6 are taken to be a trivial extension of the transverse space, with the 
6-dimensional plane wave metric being simply ds^ = 2dudV + Cabiu)dx°'dx'' with a,b = 1, . . . , 4. 
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where in MS, Z = 1 + ^(^^(^ze ^os^vr — 7). Note that this introduces a scheme 
dependence into though of course physical resuhs must be independent of this 

choice. We therefore find the following expression for A{s): 

where the renormalized kernel in MS is: 

.2 ri 



(5.36) 



The physical mass renormalization condition we are using ensures that the renor- 
malized kernel satisfies the condition S(||^) = 0. In the following section on curved 
spacetime, we will find it convenient to use a RG scheme where the freedom to add a 
further finite counterterm in Z is used to impose the additional condition S'(||^) = 0. 
This is achieved by the definition^^ 

.w^..w-..(i!)-(.-i!).'.(|!). 

The spectral function, which is independent of renormalization issues, is identified 
as before, and a short calculation gives: 

1. ^/ly^ 



piv) = Imsf 



2 . Am\m (^-^^^ 



12(47r)3 ' ' V z/2 

This time, because of the extra z/^ factor in (5.38) compared to (5.15), the relation to 
the kernel is via a twice-subtracted dispersion relation, 

t(^) = m + u^t'iO) - .3 r % P^^^ ^ . (5.39) 

This may again be checked explicitly by performing the integrals over v and ^ in (5.39) 
and (5.36). 12 



""^""^We assume here that the subtractions are real, which will be the case below threshold where 

< Am?. In general, the subtractions are specified as the real parts of and S^(|£j). 

-'^^ The integral over the spectral function on the rhs of (5.39) is explicitly 

-— — arctan ^^^= + - - z 
6 (47r)3 V ' ^/T^ 3 

for z > 1, where here z — 



-37- 



The first indication of problems witli Type I initial conditions is tlie remaining 
presence of tlie Z factor in tlie numerator of (5.35). Essentially, this is indicating that in 
a theory requiring field renormalization, and in the absence of a physical short-distance 
cut-off, the initial value at mq = must be divergent if the long-time evolution is to 
remain finite. The analytic structure is the same as for the Ac/)^ theory in c? = 4, and the 
inverse Laplace transform again separates, below threshold, into distinct .Apoic(M) and 
J^cntiu) contributions. For Type I initial conditions, the analysis of the pole contribution 
is very similar to the 4-dim theory. The new problem arises with the cut contribution 
where, substituting the spectral function, we find to C(e^): 

Contrast this with the corresponding expression for in c? = 4, eq.(5.23). Here, the 
V integral is log divergent for large z/^, due to the extra power dependence in piv). 
In turn, this can be traced back to the behaviour of the kernel for large s, where from 
(5.36) we see that S(s) ~ slogs rather than S(s) ~ logs in the c? = 4 theory. 

To evade these problems, which are intimately related to the divergent short dis- 
tance behaviour, we choose instead to analyse the theory using Type II or Type III 
initial conditions. This shifts the divergent physics off to the switch-on surface at 
Uq — > — OO, so that we start the initial value problem from m = with a renormalized 
field, already at least partially dressed. 

The Laplace transform of the equation of motion (3.18) for Type II initial conditions 

is: 

--2iujZ[sA{s) - A{^)) + M|zi^(s) + SB(s)i(s) -^(O)-(Sij(s) - Sb(0)) = , 

(5.41) 

Solving this, and writing in terms of renormalized quantities, we have 



yi(s) 



yi(o) 



S(0) + M' 



2 



2iujs - M2 - S(s) 



(5.42) 



It is straightforward to check directly in this expression how the subtractions (5.33), 
(5.34) implement the usual mass and field renormalizations. In particular, note how 
all factors of Z are absorbed by the subtraction (5.34). 

Once again, this has a pole at s = sq (see (5.6)) and a cut with branch point at 
2ibjs = Am? from (5.36). Note that, for ^ 0, there is no pole at s = despite the 
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1/s factor in (5.42). The inverse Laplace transform is evaluated as usual, and we find 



A{u) = j^poio(w) + A 



cut 



U 



(5.43) 



where now 



A 



and 

Acutiu) 



pole 



2-^(0) (l + ^) M 



u 



A{0)e- 



1 I S{0) 



1_ 1 E'(M1^ 



(5.44) 



du 

2m ^ 



p{u)e- 



■ V u 

' 2lj 



[^2_M2-ReS(|£)]VvrW^ 

factor, which affects 



(5.45) 



S(0) 



Note the differences from the Type I case - the extra 

the asymptotic ratio A{u)/A{0), and the extra M^/z/^ factor in the cut contribution, 
which plays the role of a subtraction in a normal dispersion relation in providing the 
necessary convergence factor for the integral over u. 

It is instructive to evaluate these expressions at C(e^). For the pole contribution, 
we have: 



where 



/ \ /I / \ _ i M u 

Apo\e[u) = A{0)e 2^ 
= A{0)e~— 

h(z) = 3z\/ z — 1 arctan ■ 



M2 2iuj \2iuj 



12 (47r) 



4m^\ 



3^ + 1 . 



(5.46) 



(5.47) 



Since h{z) < for all 2 > 1 (i.e. below threshold, < 4m^), the effect of the 
pole contribution is to give a w-independent reduction in the amplitude. Note that 
the renormalization scheme ambiguity cancels in the combination S(0) + |fj^'(||j), 
as is easily checked from (5.34); in particular, this ensures the disappearance of any 
/^-dependence in h{z). 

At C(e^), substituting the explicit form (5.38) of the spectral function, we find the 
cut contribution is 

4m2\3/2 



^cut(«)=^(0)^^M^ 



V 



2^ (z/2 - M-^Y 



1 



%v u 

e 



(5.48) 
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For u = 0, the integral over u can be performed analytically and, together with (5.44), 
confirms the sum rule ^poie(O) + ^cut(O) = ^(0) required for consistency. 



the amplitude oscillates on a scale uq = ^^n^'^M^ ' while the power law decay is on a 
scale md = tt^. Although the field at the initial value surface is renormalized, rather 



The physical picture emerging from (5.46) and (5.48) closely resembles that already 
described for the A(f)'^ theory in c? = 4. Evaluating the cut contribution for large u, we 
find 

Aut(n) ~ (^^j e-— , (5.49) 

similar to (5.24) but with a faster power law decay. The overall behaviour is similar to 
Fig. 3 in the below threshold case. Once again, factoring out the overall phase e~^^, 

im'^-AP ' 

A/2- 

than bare as with Type I initial conditions, once the source holding the amplitude fixed 
at 7l(0) for all u < is removed, the field relaxes in real time, tending asymptotically to 
the constant value given by Apoie{u). This further, finite, dressing of the field, as well 
as the transient behaviour for small u is of course consistent with the general unitarity 
constraints encoded in the optical theorem. 

It is also instructive to analyse this theory with the Type III initial conditions 
introduced in section 3.1. The Laplace transform of the equation of motion in this case 
is (compare (5.41)): 

-2iujZ{sA{s) - ^(0)) + M|Z.A(s) + Sb(s).A(s) 



where the inclusion of the phase with the renormalized mass makes a crucial modifica- 
tion to the final term. The solution is 



Ms) ^ ^ 



1 + 



2iujs - M2 - Us) 



(5.51) 



Again, we see how all the renormalization counterterms are absorbed into the subtrac- 
tions in the renormalized kernel. 

However, for flat spacetime - though not in curved spacetime as we shall shortly 
see - there is a further simplification. Here, it follows immediately from the mass and 
field renormalization conditions (5.33) and (5.34) that S(|ij) = 0. So then, we simply 
have 

Ms) = (5.52) 

2iu! 
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and only the simple pole remains. The inverse Laplace transform gives 



A{u) = ^(0)6-"^ . (5.53) 

The physical explanation is simple. Since for these Type III conditions, the field at the 
initial value surface has been prepared in a fully dressed, renormalized state, it simply 
continues with that evolution for m > 0. This is, however, special to flat spacetime. As 
we see in the next section, the subsequent evolution in curved spacetime can be highly 
non-trivial. 



5.4 Quantum electrodynamics 

Our flnal example is quantum electrodynamics which, as well as requiring a UV diver- 
gent fleld renormalization, is not asymptotically free. As we shall see, this introduces 
further problems with Type I initial conditions. We quote results for QED with both 
scalar and spinor "electrons', with non-zero mass m. 

The source-free equation of motion for QED in curved spacetime takes the form 




(5.54) 



In a plane-wave background, the solutions of the classical Maxwell equations are 

^^^^ = 5^.aE\{u)%{x) . (5.55) 

Here, Eia{u) is the zweibein for the transverse metric Cab{u) deflned in (2.5). The 
index, i = 1,2 (which are associated to the transverse Brinkmann coordinates) labels 
the two physical polarization states. To solve the initial-value problem, we then make 
a similar ansatz to the scalar A(f)'^ theory (compare (3.11)): 

A^(x)=A,(v^<l>l^l. (5.56) 



Specialising to flat spacetime for the remainder of this section, where the polariza- 
tion dependence is trivial and Aij{u) = A{u)Sij, the full equation of motion reduces to 
the following equation for a single complex amplitude A{u): 

-2iujA{u)+ du'J:{u,u')A{u') = , (5.57) 
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subject to the various initial conditions considered above. The Laplace transform of 
the kernel S(s) is determined from the usual momentum-space vacuum polarization 
tensor n(p^) as in (5.10) with the substitution — )■ —2ius: 

S(s) = Il{-2ius) , (5.58) 

where for scalar QED, in dimensional regularisation, 



^ (1-20' log (- 



_67r(d-4) 47r 7o ^ 47r/i2e 



(5.59) 



Note again that we are free to use the usual Feynman vacuum polarization here since 
it gives the same result in (5.57) as the Schwinger-Keldysh form. 

Gauge invariance ensures there is no mass renormalization and the field renormal- 
ization in MS is implemented by the subtraction 

t{s) = Sb(s) - {Z -l)2iujs 

= 2^u;s-J^ dai-20 log ( j . 



13 



Note that S(0) = 0. The spectral function is readily evaluated as before and we find 

''M = -?'"K5^) = m''(''-2'"'''(^-— ) ■ 

With Type I initial conditions, the solution to the initial value problem (compare 
(5.35)) is given by the inverse Laplace transform 

A{s) = ^ . (5.62) 

2ius — S(s) 

Once again, the analytic structure of the solution consists of a single particle pole, this 
time at s = 0, and the 2-particle threshold cut starting at 2ius = —Am?. In addition. 



The equivalent results for spinor QED are 



2n f n 
I](s)= 2iujs — / dfe(l-01og|- 
7^ Jo ^ 



and 

2m2\/ 4m2\i/2 



p{v) = —e{v - 2m)v' (^1 + _ ) ( 1 _ 
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though, there is the infamous Landau pole s = sl determined by the solution of the 
equation, 

l-^j\ai-2O'lo, C''^'^'f^'-^'> )=0. (5.63) 
Notice that this pole is at a non-perturbatively large Euclidean value of the momentum 

The contribution from the simple pole at s = is just a constant 

Again, the explicit presence of Z the fundamental problem of Type I conditions in a 
theory requiring a UV divergent field renormalization, although once again we see that 
to O(e^), the right-hand side of (5.64) is scheme independent. 

The cut contribution, at 0{a) is 

A,^t{u) = 2A{0) — p(z/)e-2^. (5.65) 

J2m. ^ 

With the spectral function (5.61), p(z^) ~ z/^ for large u and the cut contribution is 
logarithmically divergent. The power counting responsible for this is of course linked 
to the presence of UV divergences. 

Finally, the contribution from the Landau pole is the rapidly oscillating function 

^lp(m) = ^ , (5.66) 
a 

since is purely imaginary. This contribution is non-perturbatively large and con- 
taminates the solution for A{u) at large u. This reflects the fact that not only does 
QED require a UV divergent field renormalization but it is not asymptotically free, so 
the short- distance physics is non-perturbative. 

These problems are specific to Type I initial conditions and arise because the 
interaction is being turned on at the initial value surface, uq = 0.^^ Instead, we can 
analyse the theory with initial conditions with uq — — oo. Since gauge invariance 



^'^ Another way to see how problems arise, not specific to flat spacetime, is to recall from (3.36) that 

Q^^\u)^-— ( du" [ du' T.{u" ,u') . 

2w Juo Juo 
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ensures the photon in QED is massless, in this theory there is no distinction between 
Type II and Type III initial conditions. From (5.42), we immediately have the solution 

^ ^ s \^ 2za;s-S(s) J ^ ' 

However, in fiat spacetime (though not in curved spacetime) S(0) = and so the 
solution of the initial-value problem is trivial: A{s) only has the simple pole at s = 
and A{u) = A{0) for all u. Even the Landau pole is absent. The photon field is set 
up at the initial value surface u = already in a renormalized, fully-dressed state. No 
further dressing can take place in fiat spacetime, so the evolution for m > is trivial. 
However, this initial condition will be particularly appropriate in the following section 
when we analyse QFTs in curved spacetime, since it allows us to distinguish clearly the 
effects of curvature on dressing from the short- distance transient phenomena present 
even in fiat spacetime. 



6 Field Propagation in Symmetric Plane Waves 

We now return to curved spacetime and in this section consider field evolution in 
a special class of plane- wave spacetimes where the metric is w-translation invariant. 
These are the symmetric plane waves, or Cahen-Wallach spaces. 

The M-translation symmetry means that the intial value problem can be analysed 
using all the formalism of the Laplace transform method discussed in the previous 
section. We also note the relation with the general formalism for the refractive index 
and optical theorem described in sections 3 and 4 and our previous worlc [3, 4]. 

6.1 The inital value problem and renormalization in curved spacetime 

Since we are interested here in the effects of curvature on the evolution of the field 
and not on transient phenomena, we consider initial conditions with uq — )■ — oo, in 
particular Type HI. 

In QED, the imaginary part involves the integral 

which is divergent. In contrast, for the scalar A(j)^ theory in d = 4, the same integral has one less 
power in the denominator and is convergent. 
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The Laplace transform analysis of the equation of motion goes through exactly as 
before and so, for massive scalar A(f)'^ theory, we find 



A{s) 



A{0) 
g _ Ml 

2ibj '- 



1 



2ius - M2 - Sfsl 



(6.1) 



written entirely in terms of the renormalized quantities. The subtractions which imple- 
ment mass and field renormalization are the same as in the fiat spacetime theory since 
the UV divergences are curvature independent. However, this has important implica- 
tions for the final results for A{u). Implementing the physical mass renormalization and 
the scheme choice (5.37) for the field renormalization (adapted for curved spacetime) 
the required subtraction in the kernel is 



flat 



S(s) = Sb(s) -ReS^' 



M2 
2iuj 



ReS 



2iuj/ 

^flat / ■ 



flat/ 

B 



M2 
2ico 



(6.2) 



While this ensures the renormalization conditions ^ (jjZ;) ~ ^ ^ ^'(lIZj) = iii 
flat spacetime, crucially this is no longer true in curved spacetime and as a result the 
expression (6.1) for A{s) for Type III initial conditions becomes non-trivial. 



The simple pole in the inverse Laplace transform 



A{u) = A{0) 
is shifted to Sq, where now 



c+ioo 



ds 
2Tii 



2iuj 



2iuso = + El 



1 + 



■ 2iuj ' 



2ius - M2 - E(s) 



Evaluating the pole contribution, we therefore find 

A{0)e-—e~2c.^\2ii.)'' 



A 



pole 



.U 



2iu 



2iu 



+ ... 



(6.3) 



(6.4) 



(6.5) 



to 0{e^) in the pre-factor. As before, the cut contribution goes to zero for large 
u, so (6.5) gives the asymptotic solution for the amplitude A{u). This should be 
compared with (5.28) for A{u) in flat spacetime above threshold. Clearly, (6.5) opens 
up the possibility of particle decay in curved spacetime even below threshold. Also 
note the similarity with (3.34), which shows that the Laplace transform method has 
automatically performed the DRG resummation of secular terms into the exponent. 



In what follows, we use the leading order approximation 



A{u) = A{0)e- 



2uj^\ 2iu] ) ' 



(6.6) 
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to study the evolution of the field amplitude in a variety of symmetric plane-wave 
spacetimes. The identifications with the general formalism in sections 3.3 and 3.4 are 
evident. 

The kernel for A(j)'^ in = 4 is immediately read off from (3.24). The bare 
kernel is 

1 r°° dt , im-'t , , , 

- -2 (SF i I jVW)e-^e-" . (6.7) 

where we simply write A(t) for the VVM determinant A{u, u—t), which is w-independent 
The integral over t is defined by the Feynman prescription t ^ t — iO^, which ensures 
that the contour may normally be evaluated by rotating into the Euclidean region, 
t — )■ —it + 0"^. The integral is divergent at t = 0, corresponding to the usual UV 
divergence. Although this can be removed by a mass renormalization alone, we choose 
to use the scheme (6.2) and introduce a (finite) field renormalization as well. The 
renormalized kernel is then 




This expression clearly simplifies at s = so the crucial exponent in (6.6) is 
given by 



ReS(|£) = --^Re d^ -e^^^^^ "^(^-^^^ - ij (6.9) 

and 

In.S(|i) = J^/'<'« £fei(«'-?<So)'v/AW. (6,10) 

Note that there is no subtraction in the imaginary part. To derive (6.10), we need 
the property A*(— t) = A(i(:) of the VVM determinant. The location of the integration 
contour is important. According to the Feynman prescription, it lies just helow the real 
t-axis, in particular evading the pole at t = 0. 

The real part of the kernel S(||^) modifies the phase of A{u), while the imaginary 
part determines the amplitude. Note that both are explicitly w-dependent in curved 
spacetime, in contrast to the fiat space case. In our previous papers, we have studied 
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this frequency-dependence extensively for QED with massless photons, where it deter- 
mines the refractive index as described in section 3. All this analysis, including the 
all-important analyticity properties, goes through in the same way here for the massive 
A(j)'^ theory with only straightforward modifications. In particular, section 7 of ref.[3] 
gives explicit calculations for QED in symmetric plane waves which complement the 
discussion that follows and can usefully be read in parallel. 

Here, we wish to concentrate on the imaginary part of which determines 

the evolution of the amplitude itself, and in particular to investigate the effect 

of curvature on the decay thresholds. For these w-translation invariant spacetimes, the 
amplitude can only decrease for large enough u (beyond the transient region), 

in line with the general constraints of the optical theorem. This is consistent with 
positivity of the spectral density p{M; u) = — ^ Im'E(^). We have, 

\A{u)\ = |yi(0)|e-^^(*^)" , (6.11) 

with the decreasing amplitude corresponding to real A — )■ 00 decay with the well- 
defined rate ^p{M;uj). Moreover, the discussion in section 3 and 4 shows that this 
decay rate can only be non-perturbative in the curvature. However, we will find that 
for certain classes of plane-wave background, this decay can take place below the usual 
fiat-space threshold. 



6.2 Particle decay in symmetric plane-wave spacetimes 

We now study the rate of A — )■ 00 decay in symmetric plane- wave backgrounds accord- 
ing to the formula (6.11) with the spectral function p{M;u) given by 

Here, we have introduced the notation z = ^(^^f^ ~ '^^)' generalising that in section 
3.4. 

As a preliminary check, we recover the result already obtained in the fiat space 
limit, where the VVM determinant A(t) = 1. For z > 0, the contour can be closed 
in the lower-half complex t-plane. Since there are no singularities there, the integral 
vanishes and we simply find p{M) = 0, corresponding to the below threshold case where 
there is no decay. For z < 0, the contour is closed in the upper-half plane and picks up 
the pole at t = 0, giving 



p{M) = 




(6.13) 
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where ^± are the upper and lower solutions of the quadratic .^(1 — ^) — /M"^ = 0, 
and we recover (5.29) in the above threshold case. 



We now consider the three classes of symmetric plane wave in turn: 



(i) Conformally flat symmetric plane wave: 

The VVM determinant for a conformally fiat symmetric plane wave, for which 
0"! = a2 = a (with a real), is 



A(m, u) 



a[u — u 



sm crm — W 



(6.14) 



where u — u' = t. Inserting into (6.7) and rotating the contour t 

.2 ^1 



-it, we have 



a e 



2 (47r) 



-zt 



dt 



sinh(crt) 



(6.15) 



using v"^ = 2ius in the definition of z. 



Using the renormalization prescription (6.2), we can evaluate the renormalized 
kernel exactly and find^^ 



Us] 



1 



2 {Any 



m 



1 is 

2 ~ 2^ ^ \auj^{l-i) 



Re log 



m 



j +Re 



M2 



(6.16) 



The final term is just the optional finite field renormalization factor in (6.2). We 
can now see explicitly how the analytic structure of the kernel is changed in curved 
spacetime. Since the di-gamma function ipi^x) has simple poles at x = 0, —1, —2, . . . 
with residue -1, it follows that S(s) now has an infinite series of branch points at 



2im 



2(2p-l)a, p=l,2, 



(6.17) 



^^In the notation of section 3 and ref.[3], (6.16) involves 

-^(^) = ^(i + ^)-Relog^ , 

where z = ^(|(f3|y — M^). This should be compared directly with eq.(7.6) of ref.[3]. This reference 
also illustrates the w-dependence of the QED analogues of (6.16) for conformally flat, Ricci flat and 
general plane-wave spacetimes. 
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Consequently, the 2-particle branch cut of the flat space case has become an infinite 
sequence of branch cuts with branch points down the negative imaginary axis. 



It is useful to understand how this arises by considering a less direct method of 
evaluation. If we expand the denominator in (6.15) in powers of e~"^ and then perform 
the t integral, using a small t cut-off 5, we find 

and (6.16) can be recovered from the summation after renormalization. In this method, 
we see clearly the origin of the poles, which become the branch points (6.17) after the 
1^ integration. 

The analytic structure of S(||^) therefore comprises a sequence of branch cuts with 
branch points at 

= Am'(l + {2p~l)^^ p=l,2,.... (6.19) 

These new curvature-dependent branch points therefore depend on the parameter 
iij\/^/m? which, as we have frequently encountered, is a characteristic scale for non- 
perturbative phenomena for fields propagating in curved spacetime. As we take the 
flat space limit, the branch points all converge on the single threshold branch point at 
= 4m2. 

The spectral density is found from the imaginary part of S(s) by integrating around 
the simple poles of the digamma function. This gives 

^ ' n=0 

where ^ {n) is the root of the equation 

1 z/^ , , 
n = — + -, 6.21 

2 Auja Aujai{l - 

with < i{n) < i, and Uq is the smallest integer 

no e Z , no > -- + — . (6.22) 

2 4:UJa 

One can check that in the fiat space limit cr — ?■ 0, the sum over n becomes a continuum 
and recovers the result (5.15). 
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The long-distance behaviour of the initial value problem then follows from (6.11). 
The solution for \A{u) \ dissipates with a decay rate T = ^p{M) which is non-vanishing 
above a threshold 

> + 2uja . (6.23) 

There is then a series of further curvature-dependent thresholds at 

> 4m'^ + 2ua{2n + 1) , neZ>0, (6.24) 

at which the decay rate F jumps discontinuously. Note that the threshold is raised 
relative to the fiat space value and that in this case there is no below-threshold decay. 
We can therefore conclude that in the conformally-fiat plane wave backgrounds, the 
curvature is suppressing the A ^ (fxf) decays. 

(ii) Ricci flat symmetric plane wave: 

If we take ai = a and a2 = icr, the plane-wave is Ricci fiat and the VVM determi- 
nant is 

A(u,u')= . ^ '^>'; -. (6.26) 

sm a[u — u ) smh o"(m — u ) 

In this case the integrand in (6.5) has branch points on the imaginary t-axis as well as 
the real axis. 

In this case, the integral does not have a simple explicit solution as in (6.16). 
However, we can readily find the analytic structure of S(s) using the method above by 
expanding the integrand in a double sum coming from the sin at and sinh at and then 
perform the t integral: 

1 |-1 ^g-(5+(2p-2ig-l)<7)5 

^b{s) = -7^j-r^ / > Cpg-—j^ , (6.26) 

2 (4vr)2 Jo ^ z+{2p- 2iq - l)a ^ ' 

for some coefficients Cpq. It follows that in this case, S(s) has branch points at 

-i{2p-2iq-l)a , p,q = 1,2,.... (6.27) 



2uj 



The spectral density can be written in the form (see (6.10) 

p(z/;a;) = --ImS( — ) = --^ / / dt e"^"* ^= (6.28 
'^^ ' TT \2icoJ 4(47r)2yo Vsin at sinh at 
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The branch points on the imaginary t-axis give rise to a non- vanishing spectral density 
p(z/; u) even below the threshold u = 2m. To see this, we evaluate the t-integral 
by deforming the contour so that it wraps around the negative imaginary axis. The 
imaginary part of in the limit 4m^ — z/^ ^ aco is dominated by the contribution 

around the first banch point at t = — i^: 



(47r)2 47r J, Je 
.2 1 .1 



(6.29) 

(47r)2 2^2^ 7o 

In the same limit the ^ integral is then dominated by the saddle-point at ^ = |: 

d^e 2<T^e(i-€) ~ . / -e . (6.30) 

V ^m^ 

Consequently, in this limit we have 



.2 



As a consequence, the large u behaviour of the solution A{u) of the initial- value problem 
dissipates with a characteristic decay rate 



valid when Aim? — ^ era;. 

So for a Ricci-flat plane-wave background, the curvature induces below-threshold 
decays of A — )■ 00, with a decay rate which is non-perturbative in the curvature. 

(iii) Null energy violating symmetric plane wave: 

In this example we take ai = a2 = icr so that the VVM determinant is 

2 



A(u, u) 



a[u — u 



sinh o"(m — u') 



(6.33) 



This background violates the null-energy condition so would not be considered a valid 
solution of Einstein's equations. However, we are free to consider it as an example of 
a fixed background. It also admits a valuable check of our use of the Penrose limit to 
study field propagation in curved spacetimes. This is explained in the appendix B. 
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In this case, the integrand has simple poles on the imaginary t axis at t = —i^, 
TT G Z. In the below threshold regime, z > 0, the contribution is from the poles at 
t = —i^, n ^ Z > 0, while in the above thresold region, z < 0, the contribution is 
from the poles at t = i^, n E 'Z > 0. In both cases we can sum up the contribution 
into a common analytic function valid for either z < ot z > yielding an expression 
for the spectral density which is valid for all u: 

P^""'^ ^) = l-ihy I ^(^..A ■ (6.34) 

It is straightforward to check that the flat-space limit is correctly recovered since as 
cr — i- the integral only receives support from the region < ^ < 1 — C,o, where is 
the smallest root of the equation ^(1 — — m^/u^ = 0. 

This expression for the spectral density p(M; u) shows that the decay A — )■ 00 can 
take place even below threshold, with a rate once again non-perturbatively small in the 
curvature. 



6.3 Quantum Electrodynamics 



Now consider quantum electrodynamics in a symmetric plane wave background. The 
bare kernel for scalar QED^® was evaluated in ref. [3], and it follows directly that the 
renormalized, Laplace transform kernel is: 



a 



-u 



TC 



dt 



im t 



If. 2^5(1-0 



A,,-(t)v/A(t)f 



-St 



-I- St 



(6.35) 



using the renormalization condition (6.2). This ensures S(0) = S'(0) = in fiat 
spacetime, though not in curved spacetime, where 



a 



Re Sij(O) = —uj 

TT 



t2 



im?t 



^a(l-0 / 5^e"2-C(i-C)[A,,(t)v^-5,,] 



(6.36) 



and 

ImS,,(0) = --a;y^ da(l-0 j ^ e"2^?(W) A,,(t)v^ , (6.37) 

^''The equivalent expression for spinor QED is given in eq.(5.28) of ref.[4]. Note that in refs.[3, 4] 
there is an overall sign error in the refractive index for scalar QED, arising from omitting the relative 
minus sign between the scalar and spinor loop in the vacuum polarization (corrected in arXiv versions) . 
In particular, as we have proved here, for the symmetric plane wave examples, Imn(a;) is always 
positive, whether for scalar or spinor QED. 
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With Type III initial conditions (equivalent to Type II for massless photons), the 
field A{u) is given by the inverse Laplace transform as 



Aiu)=AiO) --— ! + - ^-^ , (6.38) 




where we have suppressed the matrix indices. This is regular at s = 0, the curvature 
shifting the simple pole to 

.0 = T^S(O) . (6.39) 
The contribution from the pole is therefore 

^poie(«) = yi(0)e2i^(°) [l - ^S'(O)] . (6.40) 
The contribution from the Landau pole is now, however, highly supressed: 

>Alp(m) ~ e--e^l^^l" , (6-41) 

so this choice of initial conditions gives a formulation of the initial-value problem that 
circumvents the problems that arise from the Landau pole and leads to a consistent 
perturbative expansion. Once again, for large u the contribution from the cut goes to 
zero, so asymptotically (6.40) gives the full result for A{u). 

Recalling the series of formulae in section 3.4, we see by comparison with the 
exponent in(6.40) that the refractive index is given to leading order by 

1 ^ 



nij[u) 



'^..-^^,,(0) . (6.42) 



That is, 



where 



a 1 

nij{u) =6ij — / rf^ ^(1 - J^iji^) 

71" Jo 



(6.43) 



poo JJ. 

7,,[z) = / - ze-^*[A,,(t)v^ - 5.,] . (6.44) 
Jo 

This reproduces the result found originally in ref. [3]. In that paper, the analytic prop- 
erties of n{u!) are explored and an extensive discussion is given (see especially section 
7 of ref. [3]) of the frequency dependence of the refractive index in symmetric plane 
wave backgrounds, showing how conventional dispersion relations are violated in curved 
spacetime while causality is maintained. Just as for the scalar A(f)'^ theory described in 
the last section, we find examples of backgrounds, notably the Ricci fiat plane waves, 
where Imn(co') = ^p{0) is non- vanishing, showing that the curvature induces the 
(necessarily below threshold) decay of the photon into electron-positron pairs. 
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7 Homogeneous Plane Waves and Singularities 



Finally, we consider the initial value problem in a class of backgrounds without u- 
translation symmetry. In this case, lm.n{u\u:) can be negative and it is interesting to 
see explicitly how this is reconciled with the optical theorem. In particular, we shall 
consider singular homogeneous plane waves, which we have already studied in detail in 
ref. [4]. These arise as Penrose limits in the near singularity region of black holes and 
in cosmological FRW spacetimes with an initial singularity. 

The profile function for a singular homogeneous plane wave (see section 2) is 
hij{u) = ^—^^5ij for some constants a^. The characterise the nature of the sin- 
gularity and display a remarkable universality [4, 12] due to their relation with the 
Szekeres-Iyers classification of power-law singularities [34, 35]. As an example, the 
near-singularity Penrose limit of a null geodesic with non-vanishing angular momen- 
tum in the Schwarzschild metric has this form with ai = | and 0^2 = |. This implies 
hij = diag(l, —1). Note that, like the original spacetime, the Penrose limit is Ricci 
flat. In general the parameters are real and > 0, or imaginary. 

The VVM matrix can easily be extracted from the equation for the geodesies in 



Since the equation is homogeneous, the solution for the geodesic spray z'^{u) = Aij{u, u'), 
defined above (2.14), is 



(2.13): 




z' = . 



(7.1) 



Aij{u,u') = «ri(W)(^-"»)/2(M"' 



(7.2) 



which gives the VVM matrix as 



ai{u — u'){uu') 2 




(7.3) 



As noted in section 2, the singular homogeneous plane waves have an enhanced scaling 
symmetry, as a result of which the VVM matrix is a function only of the ratio r = u' /u. 
In particular, the VVM determinant is simply 



where p = 1 — ("1+°^) ^ jNjQ^g p < 1 for > 0. The leading behaviour for large and 
small r follows immediately: 

A(r) ~ y/a^ rP (r ^ oo) , A(r) ~ yala; r'^ (r 0) , (7.5) 

while A(r) = 1 in the limit r — )■ 1. 

We shall focus here on the imaginary part of the refractive index. Recall from 
sections 3.3, 3.4 that this is related to the field amplitude through the formula A{u) = 
A{uo)[l + iQ^^\u)) (before applying the DRG), where 

lmQ^^\u) = u du"lmn{u";u) . (7.6) 

Juo 

We consider both scalar A(j)'^ theory and QED which, as we see, exhibit distinct and 
interesting physical effects. 

A(j)'^ theory: 

For the purely scalar A(j)'^ theory in d = 4, 



1 



(7.7) 



where 



Im J^(m; z)=lm / -e-''W ^iu,u - t) . (7.8) 
Jo ^ 

As usual, z = 2b ( g(i^g) "~ restrict to the below-threshold case where z is 

real and positive. Note that here we have kept Uq explicit, as in (3.36). In general, the 
refractive index n{u] Uq; u) is now a function of both u and the initial value Mq, though 
to simplify notation we suppress the Uq dependence and just continue to write n{u; u) 
as before. 

The analysis of cosmological FRW spacetimes requires us to start from an initial 
value surface u = Uq with uq finite and positive, rather than letting mq — )■ — cxd as we 
can do in the black hole case. This raises some important subtleties, especially with 
regard to renormalization, which we consider first. 

It follows from the description of renormalization in the context of the Laplace 
transform method in section 5 that the vacuum polarization Il{u, u'\ u, p), incorporating 
the subtraction corresponding to mass renormalization, is: 

II{u,u]uj,p) = IIb{u,u';uj,p) + 5M'^5{u,u) . (7.9) 
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From the definition = — ReS^^*(||^), together with the exphcit expression (6.7) 
for Eb(s), we therefore have 



Il{u, u'; uj,p) 



1 e' 



2 (47r) 



u — w 



-iz{u-u') ^^(^ 



The refractive index itself is 



n{u; u) = 1 



du' n(M, u'] OJ,p) 



(7.10) 



(7.11) 



Mo 



and so, changing variable to t = m — m' in the first term, we find 



n(u; to) = 1 + 



(47r)2 4w2 



di 



U — UQ 



dt 



-e-^VA(u,M-t) -Re / - 



L^o 



dt 



-izt 



(7.12) 



The UV divergence at t = cancels between the two terms in (7.12). The subtrac- 
tion, by definition, is real and so the expression (7.7) for \m.n{u]ijj) is unaffected by 
renormalization. The crucial point, though, is that the upper limits on the t-integrals 
of the bare and subtraction terms in (7.13) are not the same for finite Uq. In the fiat 
spacetime limit and with mq — > — oo, we confirm as expected that n{u\ u) = 1. However, 
even in fiat spacetime, this difference in limits implies that n{u;uj) has a non-trivial, 
M-dependent behaviour when we start from an initial value surface at finite uq. In turn, 
this implies a non- vanishing Q^^\u) and w-dependent evolution of the amplitude A{u). 

Of course, this is precisely the transient behaviour we have already studied using the 
Laplace transform method in section 5. It is interesting, however, to see this behaviour 
reproduced by the general expression (7.12) in the case of fiat spacetime, and in fact 
very similar results follow in the FRW cosmology considered below. 



We can evaluate the integral explicitly to give 

„2 1 fl 



17 



lmn(u: u) 



1 



(47r)2 4u;2 



d^ 



10 



TT 

L2 



Si{u — uq] 



(7.13) 



Note that the integrals require care near the pole at i = 0. The correct prescription matches that 
explained following (6.7), with the t contour lying below the axis and picking up a contribution from 
the pole. For example, 



Im 



-e = - — 



— sin(2t) 
1 t 
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\Aiu)/A{uo)\ 




Figure 4. The left-hand diagram shows lmn{u;u}) as a function of u in flat spacetime. 
Note that lmn{u;uj) is non-zero at the initial value surface uq and can take negative val- 
ues. The right-hand diagram shows the evolution of the field amplitude showing the 
characteristic transient behaviour as the bare field becomes dressed in real time. 



introducing the notation u = zu. Si is the sine integral. Similarly, 



Im Q{u) 



2 (47r) 



- M2 



(7.14) 



- ^o)(| - Si(n - no)) 



cos M 



{U - Mo) + 1 



Plots of Im n{u\ u) and the amplitude \A{u) \ = \A{uo)\{l - Im Q^^\u)) are shown 
in Fig. 4. These display all the features previously described in section 5 - an oscillatory 
transient behaviour where |^(m)| falls as the bare field becomes dressed, before settling 
to a fixed value in the asymptotic large-w region where lmn{u;uj) — ^ 0. Note that 
Imn(M;a;) is locally negative in the transient phase, corresponding to the oscillatory 
behaviour of 1^1(^)1 , while its integral Im Q^^\u) remains positive, in accordance with 
the optical theorem. 

Quantum Electrodynamics: 

The analysis for QED follows the same lines, but there are important differences, 
both technical and in the resulting physics. In scalar QED, there are two contributions 
to the vacuum polarization, corresponding to the Feynman diagrams in Fig. 5. 

For the Schwinger-Keldysh (retarded) vacuum polarization, we have previously 
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Figure 5. The two Feynman diagrams contributing to the vacuum polarization in scalar 
QED. 



shown [3] that 

a r p-iz(u-u') 

TT Jo 



[u — uY 

+ i6ij / — e"*"* 6{u,u') 
Jo * 



(7.15) 



where the term proportional to 6{u, u') comes from the first Feynman diagram in 

2 

Fig. 5, and as usual z = 2ujI{i-^) ■ "^^^ ^^'^ transverse polarizations are labelled by 
the Brinkmann coordinates i,j = 1, 2 and Uij can be taken as diagonal with no loss of 
generality. 

Given the relation (3.38) of the refractive index to the vacuum polarization, 

1 r ^ ~ 

nij{u;u) = 6ij - — du IlsK,ijiu,u;u,p) , (7.16) 

we therefore have 

n,j{u;uj) = Sij — / rf^^(l -0 J^ij{u;z) , (7.17) 

TT 2u Jo 

with 

ru-uo Jj. roo i. 

T,j[u-z) = / - ze-''''Aij{u,u- t),/A{u,u- t) - % / - ze"*^* . (7.18) 
Jo t Jo * 

Notice that the limits on the t-integrals in the two terms are different in general. Only 
when we take Uq — )■ — oo (as in our previous papers) do they both become the same, as 
the dependence on the initial value surface disappears. 

It is clear that the second contribution acts very much like the mass counterterm 
subtraction for A(j)'^ theory in (7.12), but with important differences. Although it 
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removes the potential singularity at t = 0, it is not a counterterm and is not simply 
the real part as in (7.12). It is essential in QED to maintain gauge invariance, and it 
is gauge invariance which keeps the photon massless, with no mass renormalization. 
Also note the t"^ power in the integrals, which arises through power counting and is 
ultimately related to the need for field (wave-function) renormalization in QED. These 
differences make a subtle, but vitally important, difference in the behaviour of the field 
amplitude in QED, even in flat spacetime. 

Consider first QED in fiat spacetime. Combining the two terms in (7.18), we have 
J^ij{u; z) = F{u] z)6ij where 

Hu: z)= § ^e-^ = / I .e-^'^--)* , (7.19) 

with the change of variable t = {u — Uo)t. Notice the occurrence of the vital factor 
l/{u — Uq), which appears because of the different power-counting in QED compared 
to A(f)'^ in d = 4. For the imaginary part (note that both vacuum polarization diagrams 
are contributing), 

lmj^(u;z) = —- / — cos(z(u — Uo)t) . (7.20) 

The integral can be done analytically and we find the following expression for the 
refractive index 



a m 



2 /-i 







TT 1 



lmn(u]uj) = / d^ SUu — Uq) 1 cos('U — Mq) , ('''•21) 

vr 4a;"' o 



2 U — Uq 



where u = zu. 



This is plotted in Fig. 6. Again we see the characteristic oscillations which occur 
when the initial value surface is taken at finite uq. The key point, however, is that 
for QED, \m.n{u\u:) diverges as l/{u — Uq) for u close to the initial value surface. 
For small {u — Uq), lmn{u;u)) is positive, so this is consistent with the observation 
in section 6 of a divergent initial dressing in a theory like QED which requires wave 
function renormalization. 

We can see this more explicitly by comparing the amplitude A{u) at some small 
value of u with its value A{ui) for some reference value Ui sufficiently removed from 
the initial value surface. Evaluating Im Q[u) from 

lmQ{ui) —lmQ{u) = u / du" lmn{u";u) , (7.22) 

J u 
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lmn{u] u) 



\Aiu)\ 



divergence 




u 



divergence 




u 



Uq 



Figure 6. This shows (left) lmn{u;uj) for QED in flat spacetime as a function of u. Note 
that for u ~ uq, lmn{u;uj) is positive and divergent. The right-hand diagram shows the 
equivalent behaviour for the amplitude |yi(n)|, showing the expected divergence on the initial 
value surface. 



with lmn{u;u) as in (7.21), we find 

\Aiu)\ ^ 



with 

lmQ{u) — lmQ{ui) 



\A{ui 



a 
2^ 



{lmQ{u) - lmQ{ui)) , 



(7.23) 



d^^il- 0\ (u- Mo)Si(M - no) + Ci{u - uq) 



+ cos[u-uq) - -{U-Uo) 



M — 7- Ml 



(7.24) 



This is shown in Fig. 6. For small {u — uq), the dominant term arises from the Ci 
function, which diverges logarithmically, and we find 

\Aiu)\ 



1 



a 



log(u - Uq) 



(7.25) 



\A{u,)\ 

We can exponentiate this using the dynamical renormalization group and conclude that 
for small {u — Uq), 

\A{u)\ 



\A{ui 



[u - Uo) 



\2-K 



(7.26) 



The physical picture, as inferred in section 6, is that in order to obtain a finite evolution 
of the photon field in QED, we have to start from a divergent value on the initial value 
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surface itself. This is the expected behaviour for a theory requiring field (wave-function) 
renormalization. Nevertheless, we can still study the evolution of the field amplitude in 
QED, either by simply taking Uq — t- — oo as in our previous work, or by comparing the 
ratio of |./1(m)| to its value at a fixed reference point away from the initial value surface 
where the interaction is assumed to be switched on. Of course, this is the essential idea 
behind the Type II or III initial value conditions introduced in section 6. 



7.1 Cosmological FRW spacetimes 

As a first example, consider a spatially- flat FRW spacetime. As shown in ref. [4], the 
Penrose limit for null geodesies with no angular momentum in the transverse space is 
a conformally-fiat, singular homogeneous plane wave with profile function^*^ 

Spatially flat FRW spacetimes are special in that the Penrose limit is a singular homo- 
geneous plane wave for all u, not simply in the near-singularity limit m ^ 0. 

We are considering an initial value problem where the interaction is turned on at 
some finite value Uq > 0, since clearly we have to start away from the initial singularity. 
We then study the development of the field amplitude along the null geodesic as the 
universe evolves. 

A(j)'^ theory: 

Specialising first to theory with M = 0, and making the convenient change of 
variable from t = u — u' to i=l — '^ = 1— r, we find from (7.8), 

ImTiu; z) = lm [ " $ e-^^"*F(t) , (7.28) 

Jo-ie t 



where z = „ J7i c\ and 



at(l-t) 2 



. ^ ucy.-Lj ^ ^7 29) 

^ ' 1 - (1 - t> ■ 



^^The coefficient 7 is related to tlie usual w parameter in the FRW equation of state p = wp 
by 7 = 3(i^w) ' ^^^^ T ~ 1' 5 '-"^ ^ rnatter, radiation or cosmological constant dominated 
universe respectively. The FRW scale factor is a{t) ^ V . The interesting example of the Milne universe 
considered in ref. [4] with 7 = 1 corresponds to w = — i, at the boundary between a decelerating and 
accelerating universe. The a = ai = a2 (to avoid confusion with the fine structure constant) parameter 
in the profile function is ci = | | . The null coordinate u in the Penrose limit is related to t in the 
FRW metric by u ^ . The fiat space limit is 7 = 0, corresponding to a = 1. 
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The integral requires careful treatment of the pole at t = (see footnote 17) and follows 
the prescription described following (6.7). There, in flat spacetime, the t-integral ran 
from — oo to oo with the contour running below the real axis - this ensured that below 
threshold [z > 0) the contour could be closed in the lower half-plane so that lmn{u; u) 
vanished. Here, the same ie prescription in (7.28) picks up a contribution ^iri Res(t = 0) 
from the pole at t = 0, and since F{0) = 1, we find 

lmT{u;z) = -- " ^ sm{zut)F{i) . (7.30) 
2 Jo t 



The refractive index is then 

.2 1 rl 



Imn(u; u) 



1 



(47r)2 4a;2 



dt 



Sim zu 



i)F{i) 



(7.31) 



The integrals over t and ^ can now be performed numerically for various values of the 
FRW parameter a. In fact, the function F{t) is very fiat over most of the range t = 
to 1. The refractive index therefore has the same qualitative behaviour as shown in 
Fig. 4. In particular, \m.n{u\uj) falls from an initially positive value and oscillates as 
u increases away from mq, and is locally negative. In those regions, Q^^\u) decreases 
and in turn the amplitude |^(m)| increases. In Fig. 7, we subtract the transient flat- 
spacetime result and plot the curvature-induced contribution to Im?T,(M; cu), viz. 



Im?T,(u; (jj) 



curved (47?)^ AuJ 



I di I " — sin(zut) [F{t) - l] . 
'■^ Jo Jo t 



(7.32) 



This shows clearly how the time-varying gravitational tidal forces induce both a 
local dressing and undressing of the renormalized quantum field. The amplitude A{u) 
decreases or increases as the screening due to the cloud of virtual pairs is enhanced 
or reduced. 

Quantum Electrodynamics: 

The situation is very similar for QED and the behaviour of lmn{u;u) and |^(m)| 
in the FRW background is qualitatively the same as in fiat spacetime, including the 
divergence for m ~ mq. Both polarizations have the same refractive index, since the 
homogeneous plane wave in the Penrose limit is conformally flat. 
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lmn{u; u) 

A(f)^ QED 




u 



Figure 7. The curvature contribution to Imra(n; w) plotted as a function o{ u — uq, for Acj)'^ 
and QED , for a fixed choice of uq in a FRW spacetime with a = |. Note that Imn(n; w) can 
be locally negative, corresponding to a curvature-induced "undressing" of the quantum field. 

If we split lmn{u; u) into the flat spacetime result already found and the curvature- 
dependent contribution, we find 

In'"(«;-)L„^ = -^^/'<i«(l-0 I cosM)[F(i)-l] . (7.33) 

where 

Fit) = Au{r)y^= \ ^ , (7.34) 

recalling t = 1 — r. Notice that the subtraction —1 in the integrand in (7.33) is the flat 
spacetime contribution, not the contribution of the extra vacuum polarization diagram 
for scalar QED. The i integration is well-defined here, since the properties F{0) = 1, 
F'(0) = ensure the absence of a singularity at t = 0. This result for lmn{u;uj) 
is plotted in Fig. 7. The results are very similar to the scalar ^40^ theory, with an 
enhancement for u ^ Uq for QED reflecting the usual power counting factor in 
(7.33). 

7.2 Black holes and the near-singularity limit 

A particularly interesting question is what happens to a dressed photon as it approaches 
a spacetime singularity. This can be addressed in a similar way, using the singular 
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homogeneous plane wave as the Penrose hmit of a Schwarzschild black hole in the 
near-singularity limit. For this case, we have ai = ^ and 0:2 = | and the geodesic spray 
defined above (2.14) is described by 



An{u,u') = 5(W)'/^(m1/^ - u"/'^ , 



(7.35) 



In the direction z^, the geodesies focus on the point u = which is the singularity. 
For the other direction z'^, the geodesies diverge as the singularity is approached. If 
we define the usual Schwarzschild coordinates {t, r, 6, 0) then by symmetry we can take 
the geodesies in the plane = 0. The transverse space-like coordinates are related 
to these coordinates via z^ = rsin6' and z'^ = rdr/du. Here, z^ is orthogonal to the 
plane of the orbit while z"^ lies in the plane of the orbit [4]. 

A(j)'^ theory: 

This time, therefore, we take the initial value surface uq < and study Imn(M; u) as 
u approaches zero. In this case, the appropriate integration variable ist = r — 1 = — 
and we have 



Im T{u; z) = lm ? e"^^!"!* F{i) , (7.36) 

Jo-ie t 



where 



Fit) = , I . + , . (7.37) 

^' ((l+t)"i-l)((l+t)"2-l) 

The VVM determinant function F{i) satisfies F(0) = 1, F'(0) = and F{i) ~ 
^Oi\a2 i'P/'^ at large i. Note that p for the background geometry is defined in (7.4). 

We are interested in the near-singularity behaviour of the refractive index, so now 
let Mo — —00. The integral over t can be expressed as in (7.30), but here, since the 
upper limit of integration becomes infinity, it is possible to rotate the contour in the 
complex t plane to run down the negative imaginary axis. Then, with r = it, we have 
the computationally more convenient form, 

lmj^{u-z)= — e-^l"l"G(r) , (7.38) 

where G{t) = ImF(— ir). Note that the lower limit is safe from singularities since 
G(0) = 0. At large r, G(r) = sin ^ tp/^ (l + ^(^-"^^'^("i''^^))) ^nd, crucially, 

is negative provided p > 0. For Schwarzschild, P = |- 



- 64 - 




\A{u) 



singularity 



log \u\ 



Figure 8. Plot of the imaginary part of the refractive index Im n(n;w) in Acj)'^ theory in a 
Schwarzschild background as u approaches the singularity at it = from u < 0. The right- 
hand diagram shows the corresponding result for the amplitude (plotted on an inverse 
log scale) showing the bounded behaviour following from Q{u) — >■ —constant as n — >■ 0. 



The resulting form for lmn{u;u)) is shown in Fig. 8. Note immediately that 
lmn{u;uj) < and diverges for small u. A key point then is whether the amplitude 
remains bounded, i.e. whether Im Q^^\u), given here by 



Im Q^^\u)=u) / du" lmn{u";u}) 



(7.39) 



remains finite as m — t- 0. This is shown in the second plot in Fig. 8, confirming that 
Im Q(^) (u) — 7- —constant as m — 0, ensuring that |^(m)| is bounded. 



To see this explicitly, we can rescale the integral (7.38) to get 



°° dr 
— e-^G 



z\u\ 



(7.40) 



then pick out the leading small \u\ behaviour from the leading term in the expression 
for G{t) for large r. This gives 

^ (7.41) 



Im J^(u; z) 



-y/oioi sin 



Tip 
T 



r 



p 



p 

[z\u\)'i 



In practice, since ai = | for the Schwarzschild black hole, the approach to the asymp- 
totic limit is slow. However, this limit is sufficient to show that Im Q^^'^ [u] is bounded 
and tends to a constant as !«! — )■ provided < | < 1, which is ensured by < < 1. 

The physical interpretation is as follows. With the initial value surface set at 
Uq —7- — OO, the incoming photon is already fully dressed as it approaches the small m. 
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near-singularity region. Here, lmn{u;u) is negative, corresponding to an increasing 
amplitude The photon is becoming "undressed", i.e. the gravitational tidal 

forces are stripping away its cloud of virtual (p pairs. This reduces the level of screening 
and the renormalized field amplitude increases as it reverts towards its bare state. Since 
this process cannot continue indefinitely, the increase in must be bounded, as 

we find. Once again, all this is consistent with the interpretation in sections 4 and 5 of 
the curved-spacetime generalisation of the optical theorem. 

Quantum Electrodynamics: 

Finally, we consider QED itself and look at the evolution of a dressed photon as it 
approaches a Schwarzschild singularity at u = 0. Here, once again taking uq — )■ — oo, 

lmnij{u;uj) = — / ^{1 - ^) lmj^ij{u; z) , (7.42) 

vr Zu Jq 

with 

— ie"^^* ^Aij{u, u-t) v/A(m, u-t)- 6ij^ (7.43) 

and with no loss-of-generality {u; z) = Ti{u; z)6ij. 

We use the same rescaling t = r — 1 = —t/u as above, so that 

Hu; ^) = n / ^ ^e-l"l*(F,(t) - 1) , (7.44) 
u Jo t^ 



where for the first polarization. 



FS) = , . (7.45) 

V ((! + *> -1) ((1 + «>-!) 

with gi = 2 — (^"1+"^) _ The equivalent result holds for the second polarization. It 
is straightforward to check that at small t, -Fi(O) = 1, -F'(O) = while at large t, 
Fi{i) -> v^afaa t^'^^ and F2{i) -> ^Jaial £92/2. 

These properties of Fi{t) ensure that the i integral in (7.44) is non-singular at t = 0, 
and since the upper limit in this case is infinity (because we have taken uq — > —00) we 
can again rotate the contour just as in the Acj)^ example and find 

1 rin- 

F,{u- z) = -— — e-^l"l- {F.{-it) - 1) , (7.46) 

H Jo 
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Figure 9. Plot of the functions Gi{T)/T'^ derived from the VVM matrix which determine 
the asymptotic behaviour of the refractive index in QED for the two polarizations. 

and so 

Im = / — e-^l"l^ G,(r) , (7.47) 

p| Jo 

where we define the functions Gi{T) = lm.Fi{—ir). At small r, Gi{T) = 0{t^), while 
for Schwarzschild spacetime with ai = |, a2 = |, we find for large r, 



(7.48) 



The small |m| behaviour of the integral in (7.47) depends on the form of the func- 
tions G'j(r)/r^. These are plotted in Fig. 9. Crucially, the asymptotic forms (7.48) 
show that these go to zero at large r fast enough, and we can define 



ci= I — Gi{t) = -0.155 , 

(7.49) 







00 



r2 
dr 



C2= I — G'2(r) = 0.105 . 




This is in contrast to A(f)'^ theory where the corresponding integral is divergent. As we 
now see, this produces quite different small \u\ behaviour for the refractive index. Note 
also the difference in sign for the two polarizations. This can be inferred from (7.48) 
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Im ni{u) 




singularity 



Figure 10. Plot of the imaginary part of the refractive index lmn{u; w) for the two polariza- 
tions in QED in a Schwarzschild background as u approaches the singularity at u = from 
u<0. 

and is determined by the properties of the VVM matrix, depending critically on the aj 
parameters. It now follows that the magnitude of the integrals in (7.47) will increase 
as |m| becomes smaller and the exponential decays less rapidly with r, until in the limit 
\u\ — )■ we find 

1 /"°° fir r- 
Im J-,(u; z)^-— — G,(r) = -f^ . (7.50) 

\u\ Jq \u\ 

We therefore find the following results for the small |m| behaviour of the imaginary 
part of the refractive index in QED (see Fig. 10): 



lmni{u;uj) = < ^ ' 

OTT 2UJ \u\ 

a 1 C2 
lmn2{u;uj) = ^T^-pT > ^ • 

DVT 2u \u\ 



(7.51) 



Imn(M;ci;) diverges as in the near singularity limit for both polarizations, but with 
opposite signs: polarization 2 is being dressed while polarization 1 is being undressed 
by the gravitational tidal forces. 

The most illuminating way to see the effect on the amplitudes is to compare A{u) 
with the amplitude A{ui) for some fixed ui sufficiently far into the near singularity 
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region that the above approximations hold, then let \u\ — )■ 0. From 

\Mu)\ 
\Mui)\ 

we find 



1 — uj du" lmni{u";u) , (7.52) 



+ — — Cjlog— , (7.53) 



\Ai{ui)\ 127r ^ui 

for the two polarizations. These are both logarithmically divergent, but we can use the 
dynamical renormalization group to exponentiate to give 

\Ai{u)\ _ ( U\Tt^^ ^^^^^ 



\Ai{ui)\ \ui 

Recalling (7.49) for the coefficients q, we finally find the small \u\ behaviour: 



\A2{u)\ ~ \u 



0.105-r^ 



127V 

(7.55) 



so in fact only |^i(m)| diverges, while for the second polarization |yi2('u)| — )■ as \u\ — )■ 0. 
The diverging case is for the polarization along the direction for which the geodesies 
focus. 

We see therefore that the evolution of the field amplitude as the singularity is 
approached is determined by the properties of the VVM matrix, which is encoded in 
the signs of the coefficients Cj and differs for the two polarizations. For polarization 
2, which correspond to the direction along which the geodesies diverge and for 
which lmn{u;uj) is positive, the photon becomes increasingly dressed - the screening 
increases and the amplitude falls to zero. For polarization 1, however, for which the 
geodesies focus, the gravitational tidal forces near the singularity strip away the vacuum 
polarization cloud of virtual pairs and the photon becomes undressed, the screening 
is reduced and the amplitude increases. In the purely scalar Acj)'^ theory, this process 
is bounded and the amplitude rises to a limiting value. In QED, however, we see the 
consequence of the theory requiring wave function renormalization. As the singularity 
is approached, the renormalized, dressed photon is being restored in real time to its 
bare state, but since for QED the ratio of bare to renormalized fields is UV divergent, 
in this case the amplitude rises without bound as \u\ — )■ for one of the polarization 
states. 
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8 Conclusions 



This work was motivated by the apparent paradox that, when vacuum polarization is 
taken into account, photons propagating in curved spacetime can experience a refrac- 
tive index with a negative imaginary part, corresponding to an amphfication of the 
amphtude. This is in direct contradiction to the conventional flat-spacetime optical 
theorem, which relates lmn{uj) to a decay rate and is necessarily positive or zero. 

In this paper, we have addressed this issue from the point of view of an initial value 
problem, tracing the evolution of a renormalized (dressed) quantum field as it propa- 
gates through curved spacetime. Different initial conditions were used, corresponding 
to an initially bare or dressed field, and renormalization issues carefully addressed. The 
simplifications following from being able to replace the background spacetime by its 
Penrose plane-wave limit were fully cxloited, allowing in many cases a very detailed 
analysis of the field evolution using a Laplace transform solution of the one-loop equa- 
tion of motion. 

The physical picture that emerges is that in a curved spacetime background, gravi- 
tational tidal forces act on the virtual e'^e~ cloud which dresses a renormalized photon 
field in such a way as to increase or decrease the degree of dressing. In particular, grav- 
ity can "undress" a renormalized photon, returning the field towards its bare state. As 
the degree of dressing is reduced, so is the level of screening and in consequence the 
amplitude increases. This effect can therefore produce a negative lmn{u;u), though 
the effect is constrained to be local, and bounded (with the exceptions above of theories 
with UV-divergent wave function renormalization taken at the limit of the intial-value 
surface or a spacetime singularity). Indeed, the propagation of a renormalized pho- 
ton through a time-dependent curved background resembles in many ways the initial 
transient phase of propagation in fiat spacetime when an interaction is abruptly turned 
on. 

In addition, in studying scalar A(f)'^ theories, as well as QED, we found many 
examples of below-threshold decays which would be kinematically forbidden in fiat 
spacetime. These decay rates are non-perturbative in the curvature. In particular, 
the analytic structure and nature of the thresholds was studied in detail in the class 
of symmetric plane waves, where a rich pattern of curvature-dependent thresholds 
emerged. It is probable that these can be interpreted in terms of group representations, 
since this class of homogeneous plane waves admits an enhanced symmetry described 
by a Heisenberg algebra. 
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The central result of this paper, which provides the formal framework for all of 
these phenomena, is the formulation of a generalised optical theorem valid in curved 
spacetime. In a plane-wave spacetime, this has the form: 

PA-^4,<f){u) = Auj I du" lmn{u";u) = / du" / du ImllsKiu" , u' ; u , p) , (8.1) 

and relates the total A — )■ 00 decay probability to the integral of the imaginary part 
of the vacuum polarization. The key point is that while this allows lm.n{u\uj) < 
locally, when integrated over the entire trajectory from the surface u = Uq dX which 
the interaction is turned on, JJ*^ du"lT[in{u"] uo) must be positive to ensure consistency 
with unitarity. 

This confirms the theme of our earlier work that, despite the many novel and 
unexpected effects due to vacuum polarization, quantum field theories in curved space- 
time indeed respect the fundamental principles of causality and unitarity. However, in 
so doing, many of the basic and widely-assumed theorems of QFT in flat spacetime, 
notably dispersion relations, the analytic structure of Green functions and scattering 
amplitudes, as well as the optical theorem, must be reformulated in curved spacetime. 
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Appendix A: Weak Curvature Expansions 

In this appendix, we quote the results for the weak coupling expansions of the 
refractive index for scalar and spinor QED generalizing the case of massless scalar 
fields (3.43) 

The case for scalar QED follows from (7.17): 

riijiu^uj) = ^(-R««(^)5ji + 2Riuju) 



iau 



16807rm4 



(u)) + 



(A.2) 
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Here, Ruu = Riuiu + R2u2u is a component of the Ricci tensor. The resuh for (spinor) 
QED can be extracted from [4] 



iau ■ ■ i-^-'^) 
- X2607rm^ {25Ruu{u)6ij - QRiujuiu)) H , 

The first correction here is the original Drummond-Hathrell result [10]. The higher 
corrects are then non-linear in the curvature. 

These expansions allow us to extract the leading order behaviour of the amplitude 
generalizing the expression for the scalar field in (3.44). For scalar QED 



■Ai(Mi) = Aj{u2) exp 



and (spinor) QED 



Ai{ui) = Aj{u2) exp 



iuju 



leSOvrm^ 

(A.4) 



^ {25R 

uu 

(Ml )(5jj QRiuju 

(ui) - 25Ruu{u2)Sij + QRiuju{u2)) 

(A.5) " 
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Appendix B: Penrose limit and de Sitter space 



A key element of our analysis of field propagation in curved spacetimes has been 
the simplification brought about by the Penrose limit. It is remarkable that a highly 
non-trivial validation of this method can be found using the symmetric plane waves of 
the third type considered in section 6. 

To see this, consider the spacetime dSs x M, with metric 

ds^ = -de + cosh2(at) [de^ + sin^ 6 dcj)^) + dz^ , (B.l) 

The three-dimensional de Sitter space metric is given in global coordinates where {6, </>) 
are the coordinates on S*^. A congruence of null geodesies consists of t = z with 6 = Oq 
and = 00 fixed. The associated null coordinates are 

u=^iz + t), v = l=(z-t). (B.2) 

The Penrose limit is now trivial to take: we identify = 6 — 6q and = sin~^(^o)(0~ 
00 ) and expand in powers of v and x° keeping terms of order 2 with v having weight 2 
and weight 1: 



ds' 



= 2du dv + cosh^(crM)da;" dx" . (B.3) 

Penrose limit 



with cr = (\/2a)^^. In this case, the Penrose limit is the symmetric plane wave of the 
"wrong sign' kind with ai = cr2 = icr- This geometry violates the null energy condition 
since in Brinkmann coordinates = — 2(t^ < 0, but we can nevertheless consider 
this, as in section 6, valid fixed background. 

The remarkable feature that allows the test of the Penrose limit method is that 
the Green functions and spectral density are known exactly in this case for the original 
spacetime dS^ x M, as well as its symmetric plane wave Penrose limit, due to the 
high degree of symmetry of the de Sitter space. In particular, the Green functions 
admit a Kallen-Lehmann representations from which the spectral density p(M) can be 
extracted. In the limit M — )■ 0, we will show how this exact result reproduces the 
spectral density (6.34) already found in the plane wave background. 

Using the notation of [38, 39], the Kallen-Lehmann representation is 

POO 

G^{iy; X, x)G+{p; X, x) = / dn^ py^^{K)Gj^{K]x,x') . (B.4) 
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Here, G+{iy; x, x') is the Wightman function and for a (minimally-coupled) field of mass 
m 

.^^^-1. (B.5) 

where x are coordinates on dS^. The Kallen-Lehmann weight-function for three- 
dimensional de Sitter space is explcitly^^ 

sinh^(7rfi;) 

^"'^ ~ 2%^ cosh "("+;+^) cosh cosh llMp^ cosh ■ ^^-^^ 

The Green functions for a field of mass m on the product space dSs x M follows 
in a simply way by noticing that momentum along the z direction acts as an effective 
contribution to the three-dimensional de Sitter mass. Splitting the coordinates as 
X = {z, x), we have 

^r\x,x') = I |^e^^(— ')G+(v%^-l;x,£0 . (B.7) 



GT'ix.x', . ^ 

ZTT 

Taking the Fourier transform along the z direction, we have 

(p;x,x') = G+(v/^-l;x,x') . (B.8) 

The vacuum polarization for the scalar A(f)'^ theory is 

Il{x,x')=ie^G^+\x,x'y . (B.9) 
Using the Kallen-Lehmann representation for dS^ in (B.4), we therefore find 



oo 



dq 

dM^ p{M)G^f\p;x,x') , 



dze-'^^'-''^j:{x,x') = e^ / ^Gf)(g;x,x')G'f^(p-g;5,x) 



(B.IO) 







PiM) = ^ I 7^P.,.(«:) , (B.ll) 



where the weight function is 

f dq 

with 

.^ = ^^-1, ^.^ = "^' + (^^^'-1, .^ = ^^-1. (B.12) 
2a2 ' ' 2(t" 
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This follows directly from equations (4) and (38) of 
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In particular, p(M) is precisely the spectral density for the vacuum polarization on 
dS^ X M and the decay rate for A — )■ 00, for massless A, is given by ^p{0). 

Now we carefully take the "geometric optics" and "weak curvature" limits described 
in section 2.1, which led to the use of the Penrose limit. First of all, we identify the 
external momentum as p = where u is the light-cone momentum. In the geometric 
optics limit, due to the form of the function pi,^^(/t), the integrand in (B.ll) only has 
support in the neighbourhood of g = |, i.e. where the external momentum is shared 
equally by the two particles in the final state. Consequently the arguments of all the 
hyperbolic functions in (B.6), except the one involving K — u — fi, are large and therefore 
may be replaced by exponentials. This means we can approximate 

In particular, we can parameterize Q = ^ and restrict the C, integral to the interval 
[0, 1] and then expand for large u 

2 

m 

It is important for the validity of these expansions that the support of the integrand 
lies away from the points = 0, 1. Finally, putting all this together gives the spectral 
density at M = 0: 

e2 .1 



l^^J ^0 1-1- g2^<Te(i-o 



This is precisely equal to (6.34) for M = 0. This confirms that the results we 
have obtained directly from the Penrose limit spacetimes are indeed the correct ap- 
proximations to those of the full background spacetime when we impose the physically- 
motivated conditions described in section 2. 
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